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ABSTRACT 


Dynamical systems for the nuclear fission reactors based on point 
reactor kinetics with reactivity feedback effects have been studied using modern 
concepts from nonlinear dynamics. 

The governing equations were suitably transformed so as to minimize 
the number of parameters. The possibility of the existence of limit cycles 
for the various models was then investigated by analytical methods as well 
as numerical experiments. Limit cycles were ruled out in one-temperature 
models in the neighborhood of the operating point. In the two-temperature 
model without delayed neutrons, a limit cycle was predicted using the method 
of Hopf bifurcation. This was confirmed by numerical experiments, thus leading 
to the conclusion that the predicted cycle was stable. The necessary condition 
for Hopf bifurcation was also established in the two-temperature model with 
delayed neutrons. 

Some additional investigations on the possibility of chaos were 

; /' ■/ '* . . ; : ■ 

also done by computing the largest Lyapunov exponent. The result was in the 
negative for the particular parameter values suitably chosen after the Hopf 
bifurcation in the two-temperature model without delayed neutrons. 
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CHAPTER 1 


INTRODUCTION 


To begin with, in this chapter we introduce the subject, followed 
by a review of literature and an outline of the present work. 

1.1. Prologue 

Having created a giant in the form of nuclear fission reactor, 
man has to make sure it serves him as an obedient and efficient servant, 
without becoming a threat to its own master. Safe design and control of nuclear 
reactors necessitates an understanding of the transient behavior of reactor 
state variables, particularly the neutron population and fuel temperature, 
with respect to various parameters. This is the subject of study of nuclear 
reactor dynamics. 

By mathematical modeling of a nuclear fission reactor, we get 
some equations which form a nonlinear dynamical system. Linear analysis 
of such a dynamical system only tells us about the initial effect of small 
disturbances about the operating point. In order to have a better understanding 
of the dynamics, we have to go for nonlinear analysis and bifurcation theory 
using state space concepts. Nonlinear analysis shows us more interesting orbits 
in state space than the linear analysis, which shows only stable and unstable 
fixed points, that too, only locally. Nonlinear analysis of a system may show 
periodic orbits (limit cycles), quasiperiodic orbits (tori), and aperiodic or chaotic 
orbits (strange attractors). By studying bifurcations, we may find out whether 



such orbits occur in a particular dynamical system, and if so, at what values 
of the parameters. 

If, for certain values of parameters, a disturbance from the operating 
point is going to take it to a limit cycle, ending up with ceaseless oscillations 
in neutron population, fuel temperature, etc., then such a phenomenon is certainly 
worth studying. Hence the study of limit cycles in nuclear reactor dynamics. 
A limit cycle is an isolated periodic solution of an autonomous system (i.e. 
a system without explicit dependence on time) represented in the state space 
by an isolated closed path. The neighboring paths are not closed, but spiral 
into or away from the limit cycle. The study of Hopf bifurcation provides 
us with a clue to the existence of limit cycles. In this type of bifurcation, 
a pair of complex conjugate eigenvalues of the 3acobian matrix crosses the 
imaginary axis transversely, i.e. with nonzero velocity, and makes it possible 
for a limit cycle to exist. 

A limit cycle may further bifurcate to a 2-torus, i.e., a quasi- 
periodic or almost periodic orbit with two incommensurate frequencies, by 
undergoing a secondary Hopf bifurcation. Still further bifurcations may lead 
to a 3-torus and subsequently to chaos. Moreover, a system may go to chaos 
via other routes. The study of chaotic dynamics is very complicated. In this 
work we confine our scope mainly to the study of Hopf bifurcation and limit 
cycles in fission reactor dynamics, which is of great practical importance. 
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1.2 Time Dependent Modeling of Nuclear Fission Reactors 


In nuclear reactor dynamics, a mathematical model of the nuclear 
reactor is obtained [1-4] in the form of partial differential equations describing 
the dependence of state variables, such as neutron population, reactor tempe- 
rature etc., on space and time, and parameters. There are a number of methods 
for deriving reactor kinetics equations : they range from the rather complica- 
ted to the very simple [1]. 

1.2.1. Space Dependent Dynamic Model : 

The equation for time dependent diffusion of neutrons in a reactor 


is [3] 


It- = Dv V 2 N - I vN + S 
d t a 


where, 

N (ryt)dV = number of neutrons in a volume element dV at a point jr at 
time t ; 

2 • - 

Dv V NdV = number of neutrons diffusing into dV per unit time at 

time t ; 


Z vNdV 
a 


= number of neutrons absorbed in dV per unit time at 
time t ; 


S(ryt)dV 

D 

Z 

a 

v 


number of neutrons produced in dV per unit time at 
time t. ; 

diffusion constant ; 

macroscopic neutron absorption cross-section.; 
neutron speed. 


Here, it has been assumed that the neutron current density is given by Fick's 
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law, j_ - -Dv VN, and that all the coefficients are independent of position, 
and that their numerical values represent suitable averages over the neutron 
velocity distribution. Let, 


^oo = infinite medium reproduction factor, 

= delayed neutron fraction for the i-th precursor, 

6 - 6 i> 

^ j = decay constant of the i-th precursor, 

Cj(r_,t) = volumetric density of the i-th type of precursor. 


Including sources of neutrons extraneous to the fission processes as a term 


S 0 (£?t), we get diffusion equation and the delayed neutron equation as [3]. 


3 N 
3 t 


Dv v N - 2 a vN + (l-B)k 00 Z_vN + A. C- + S (1.2) 


3C. 

i 

3 1 


3 ; k 2 vN - A. C. 
i 00 a ii 


(1.3) 


Thus we get a set of coupled partial differential equations, i.e. 
an infinite-dimensional dynamical system. To reduce it to a finite-dimensional 
system, i.e., a set of ordinary differential equations, we use the lumped para- 
meter or point reactor model. 

1.2.2. The Point Model of Reactor Kinetics : 

In this model, it is assumed that all but the fundamental spatial 
harmonics of the flux, following a flux disturbance, have been attenuated 
in a very short time (a few microseconds) compared to times of reactor kinetics 





interest [I]. That is, the neutron density in the reactor is now assumed to 
rise or fall as a whole following a reactivity change. Hence all spatial dependen- 
cies factor away, leaving time as the only remaining independent variable. 

Assuming that N and C. are separable in space and time, we subs- 
titute [3] 

N(r,t) = f(r) n(t) ; C^ryt) = g^r) c^t) (1.4) 

into Eqs. (1.2) - (1.3). Then the removal of space dependence from the resulting 

2 

equations requires that f/g^, V f/f and S Q /f be independent of position. Assuming 
g i = f; V 2 f + B 2 f = 0 

2 

(where B is the fundamental mode buckling) and letting 
q(t) = S o (r,t)/f(r), 

we get the point reactor model as 


dn 

dt 

dc. 

i 

dt 

where , 

P 

1 



I. X.C; . q 



X. 

1 


c. 


l 


reactivity, 

neutron generation, time, and 


= effective source strength. 


(1.5) 

0 . 6 ) 


At high power, we may ignore the term q in Eq. (1.5) [3]. Further, we 
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use one group representation, i.e., we assume that all the delayed neutrons 
have a single mean lifetime. Then we get the equations 

(1.7) 

0 . 8 ) 

This is the simplified point kinetic model which is a system of 
two coupled ordinary differential equations. The inclusion of feedback effects 
will make the first equation nonlinear. This will be discussed in detail in 
Chapter 2. 

1.2.3. Discrete Models : 

A further simplification in the finite dimensional continuous dynami- 
cal system of the point reactor model may be attained by discretizing the 
equations, or by starting with discrete population dynamics models. Such discrete 
models are the ones which are used in biological population studies, and frequently 
exhibit complicated or chaotic behavior. It is an open question whether or 
not discrete models are more appropriate than the continuous models in the 
study of neutron population in a nuclear reactor. We will have more to say 
about this in a later chapter. 

1.3 Literature Review 

Several papers have appeared [5-13] on the stability analysis and 
oscillatory behavior in nuclear reactor dynamics. Reator kinetics books [1-4] 


dn 

dt 


p - e 


n + X c 


dc 

dt 


n - A c 
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have also discussed these aspects of nuclear reactor dynamics. 

1.3.1. Books : 

According to Ash [1], it is of paramount importance for stability 
determination to ascertain the existence of limit cycles, however, there seems 
to be no general way to determine whether or not a limit cycle exists. Earlier 
boosk by Lewins [2] and Hetrick [3] have shown the existence of periodic orbits 
in the one-temperature feedback model by neglecting the effect of delayed 
neutrons and assuming a constant power removal during the transient. These 
periodic orbits form a nonlinear center rather than limit cycles. Akcasu et al 
[4] have shown that in some reactor models, the reactor power does not return 
to its equilibrium value following an initial perturbation, nor does it diverge 
from it at a finite or infinite time. Instead, it oscillates between two limits 
about the equilibrium power level as time unfolds. They have called such solutions 
of the kinetic equations as bounded solutions or limit cycles. 

1.3.2. Research Papers : 

Akcasu and Shotkin [5] have predicted the existence of limit cycles 
for reactors with arbitrary linear feedback and no delayed neutrons, but they 
have admitted that proving the existence of limit cycles convincingly is a 
delicate problem. They have argued that when a reactor is linearly unstable 
but bounded, closed periodic curves (limit cycles) must exist. Devooght and 
Smets [6] have considered the determination of regions of stability. Schmidt 
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of stable and unstable limit cycles in two-temperature feedback model without 
and with delayed neutrons, respectively. Shotkin et al [8] have shown that 
the presence of delayed neutrons may result in unstable limit cycles. Extending 
the approach of Akcasu and Shotkin [5], they have predicted the existence 
of an unstable limit cycle by showing that the system is locally stable but 
not stable in the large. Vreeke and Sandquist [9] have demonstrated the existence 
of stable limit cycles in two-temperature feedback model without delayed 

neutrons by phase space computations. Keener and Cohen [10] have analysed 
the same model as in Ref. [9] using a perturbation method and shown the 

existence of stable limit cycles. Poddar et al [11] have analysed limit cycles 
using perturbation techniques ignoring the presence of delayed neutrons. 

In most of the above works, the presence of delayed neutrons 

was ignored. Ashok Kumar [12] analysed the two-temperature feedback model 
with delayed neutrons. Using a perturbation approach for linearization, he 

established the conditions for oscillatory and non-oscillatory transient behavior 
as well as sustained oscillations, and verified the results of this linear analysis 
by numerical computations. Enginoi [13] has obtained a criterion for asympto- 
tic stability of a reactor in terms of the equilibrium power level and delayed 
neutron parameters. 

1.4. Outline of the Present Work 

In the literature cited, modern concepts of nonlinear dynamics 
such as Hopf bifurcation have not been used to predict the existence of limit 
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cycles and analyse their stability. The present work attempts to fill this gap. 

In Chapter 2, feedback effects are included in the equations of 
point-reactor kinetics. The dynamic equations thus obtained are transformed 
to a dimensionless form, and dynamical systems for various reactor models 
are constructed. 

The dynamical systems constructed in Chapter 2, are explored 
for the possibility of Hopf bifurcation and the existence of limit cycles in 
Chapter 3. The predictions based on the analysis are verified by numerical 
experiments. 

Some additional investigations on these systems, to see the possibility 
of chaos, are reported in Chapter 4. 

Finally, Chapter 5 concludes the thesis with a summary of the 
present work and some suggestions for further work. 


* * * 





CHAPTER 2 

CONSTRUCTION OF DYNAMICAL SYSTEMS 
FOR NUCLEAR FISSION REACTORS 


In this chapter, the point reactor model with energy balances and 
reactivity feedback is formulated. These^ equations are put in a convenient 
form by a suitable translation and normalization. Then dynamical systems 
for various models are constrcuted by making some approximations. 

2.1. The Lumped Parameter Dynamic Model 


In Section 1.2.2, we obtained Eqs. (1.7) and (1.8) for the point 
reactor model with one group of delayed neutrons. Writing reactor power, p 
instead,., of neutron density n, and using power units for c also, we can rewrite 
the equations as 


dR 

dt 


P-g 

1 


P 


+ 


X c 


( 2 . 1 ) 


dc 

dt 



( 2 . 2 ) 


To these equations, we shall now append the equations obtained 
from energy balances and include their feedback effect in the first equation. 


2.1.1. Energy Balances : 

The lumped parameter thermal model is based on two energy balances: 
one for the energy contained in all the fuel elements of the reactor core, 
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and another for the energy stored in the coolant within the core volume. 
The energy balance for the fuel elements is given by 


dT, 


dt 


h f < T f - T c> 


where, 

Cj = heat capacity of the fuel elements in the reactor 

core (3K' 1 ), 

hj = fuel to coolant heat transfer coefficient (WK - *), 

Tp T c = average temperatures of the fuel and the coolant, 

respectively. 

Dividing ail the terms by C^, we can rewrite the equation as 


dTj/dt 


= Kp - y (T f - T c ) 


( 2 . 3 ) 


where, K = 1/C^ and y = h f /Cj 


The energy balance for the coolant within the core is given by 


dT 


'c dt 


= h / T f - T c> - h c (T c - T i> 


where, 


heat capacity of the coolant within the reactor 
core (JK ')» 


T. 

1 


= coolant inlet temperature, and 


twice the product of the mass flow rate and 

- 1 . 


specific heat of the coolant (WK- ). 


Dividing all the terms by C , we can rewrite the equation as 



(2.4) 


dT c /dt = Y r(T f - T c ) - s (T c - T.) 

where, r = C^/C c , s = h c /C c , and y is as defined in Eq. (2.3). 
2.1.2. Reactivity Feedback Effects : 


The reactivity P (t) that drives a transient is the net effect of 
contributions arising from the external disturbances to which the transient 
is due and the termal hydraulic feedback effects. Therefore, it can be written 
as 


P (t) = pj(t) + p fb (t) 


where p.(t) and p^t) are the reactivities from the external disturbances 
and thermal hydraulic feedback effects, respectively. Assuming linear feedback, 
we write 


P fb (t) 


= “f (T f- 


t R V 

T f >* 


a (r - T 
c c c 


R 


where, 



fuel and coolant temperature coefficients of 
reactivity, respectively, and 


= some reference values of T^, T_ respectively. 

These reference temperatures may be chosen arbitrarily. If we choose them 
as the steady state operating temperatures T° and T° , then p. turns out 
to be zero. Therefore, for the linear feedback model, we have 



= aiT, - T° ) + a (T - T° ) 
f f f c c c 


p 


(2.5) 



For higher order feedback, we can write 


0 =c j(Tf -T°)* a c (T c -T°) 

+ F <T { - T° , T. - T° ) (2.6) 

where F (x, y) is a nonlinear function with no constant or linear terms, i.e. 

F(o,o) = o ; F (o,o) = F (o,o) = o 
x y 

2.2. Derivation of Dimensionless Equations 

Eqs. (2.1) - (2.4) appended by Eq. (2.5) represent the lumped parameter 
dynamic model with linear reactivity feedback. It has two fixed points : the shut 
down point (0,0, T., T.) and the operating point (p Q , c q , T° , T° ), where p Q and 
c q are the steady state operating values of p and c, respectively, in the final 
system (after the disturbances). Now we obtain dimensionless form of these 
equations with appropriate normalization and translation. 

2.2.1. Normalization of State Variables : 

In choosing the normalization and translation of the state variable^ 
the considerations were : 

(1) The resulting transformed state variables should be dimensionless, 

(2) The normalization should not be done with respect to parameters 
passing through zero, 

( 3 ) The number of parameters in the final equations should be minimum, 


(4) The operating point should be brought to the origin. 

(5) The shutdown point should be independent of the parameters, if 
possible. 

(6) The transformed variables and the parameters should allow some 
direct physical interpretation, as far as possible. 


Several types of normalization and translation were tried, and, 
on the basis of the above criteria, the following set was chosen : 


X 


1 


(P - P ) (C - C ) 

0 _ Y _ 2_ 

P 5 a 2 " C 

o o 



2.2.2. Time Normalization : 


(2.7) 


( 2 . 8 ) 


In deciding time normalization, the criteria were : 

(1) A dimensionless time should be obtained. 

(2) The number of parameters in the equations should be reduced to 
a minimum. 

(3) The orders of magnitude of the final parameters should be reasonable. 
(^) The time scale should not become too inflated, so as to require 

too large integration steps. 

With these considerations, several time normalizations were tried, 
and the following one, with respect to the decay constant, was chosen : 



2.2.3. The Dimensionless Equations : 


Making the aforementioned transformations in Eqs. (2.1) 
to (2.5), we get the dimensionless form of the equations as : 


— -bX j + bX£ + aX^ + qX^ + aX^X^ + qXjX^ 


X. 


X 1 -X 2 


X 3 = k(l-p) Xj - kX 3 + kpX 4 


X, 


kr 


X. 


P ^ 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 


where overhead dots denote derivatives with respect to t. This system has 
six parameters : 


b = 8/ XI , a = a f (T° - T.) / A 1, q = a c ( T c ' T i )A 1 ’ 

k = T/X = h f /C f X = KP o /X (T° - T° ), 

p > - T.) / (T° - X ) , r = C f /C c = l/KC, , 


where the various quantities appearing on the right hand sides have already 
been defined, and all the new parameters, as well as the dependent and indepen- 
dent variables are dimensionless. 

The system of Eqs. (2.10) - (2.13) has two fixed points: 
the operating point (0,0, 0,0) and the shutdown point (-1, -1, -1, -!)• This suggests 
that the values of state variables below - 1 would be unrealistic. Using the 
equations with appropriate approximations, we can construct dynamical systems 
for various types of point reactor models. 



2.3 Two-Temperature Feedback Models 


The thermal model used in the previous sections was based on 
two-temperature feedback via two energy balances. Delayed neutrons were 
included and the feedback was assumed to be linear . 

2.3.1. Model with Delayed Neutrons : 

This is the model already given by Eqs. (2.10) to (2.13). Defining 
for convenience, new parameters, 

h = kp, f = kr/p, 
we can rewrite the equations as : 


*1 

= -bXj + bX 2 + aX 3 + qX^ + aXjX 3 + qXjX^ 

(2.14) 

*2 

= x r x 2 

(2.15) 

*3 

= (k-h) Xj - kX 3 + hX^ 

(2.16) 

*4 

= fX '-fX„ ' 

3 4 

(2.17) 


This is a four-dimensional continuous dynamical system with six 
parameters and two fixed points (0,0, 0,0) and (-1, -1, -1, -1). 

2.3.2. Model without Delayed Neutrons : 

To obtain this model, we simply put b = 0 in Eq. (2.14) and forget 
about Eq. (2.15). Thus we get the model : 



17 


X 2 = (k-h) Xj - kX 2 + hX 3 (2.19) 

X 3 = fX 2 - fX 3 (2*20) 

where we have written X 2 and X 3 respectively, in place of X 3 and X^ of 
the previous system. It is understood that in the calculation of the parameter a, 

the value of 1 will be suitably changed to an effective lifetime so as to compen- 

sate., as far as possible, for ignoring delayed neutrons. 

This is a three-dimensional continuous dynamical system with five 
parameters and two fixed points (0,0,0) and (-1, -1, -1). 


2.3.3. Higher Order Feedback Models : 


So far, we have assumed the reactivity feedback to be linear. 
To include the effect of higher order feedback, we transform Eqs. (2.1) to 
(2.4) in the same manner as before, using Eq. (2.6) instead of Eq. (2.5). Thus, 
we can replace Eq. (2.14) by 

Xj = -bXj + bX 2 + aX 3 + qX^ + aXjX 3 + qXjX^ 

+ (1+Xj) tp(X 3 ,X 4 ) (2.21) 


where \p is related to the nonlinear term F in Eq. (2.6). It can be easily seen 
that 


3i}j 


i|> (0,0) = 0, 


= 0 and 


3 1 (0,0) 


7 % 


= 0 


( 0 , 0 ) 


Similarly, Eq. (2.18) (in case of absence of delayed neutrons) can 
be replaced by 


= aX 3 + qX^ + aX j X 3 + qXjX^ + ( 1+Xj) 4 1 (X 3 ,X^) (2.22) 



ID 


The models with nonlinear feedback have the same sets of parameters 
and fixed points as those with linear feedback, with the addition of a nonlinear 
function. 

2.fr. One-Temperature Feedback Models 


The one-temperature feedback model is based on the assumption 


that a, = 0 and that the heat capacity of the coolant within the core, C c , 
is negligibly small, so that, it must transfer away all the heat it receives 
from the fuel elements. Therefore, in Eq. (2.13), we take the limit as r+ « and 
obtain the relation X 3 = X^. Substituting this in Eq. (2.12), we get 

X 3 = k(l-p) Xj - k(l-p) X 3 

or X 3 = cXj - cX 3 (2.23) 


where , 


= k( 1 -p) = 


y< T 


T c °) 


Kp, 


X(T° - T.) 


X(T° - T.) 


hf(Tf ~ V 


Thus, we have eliminated one equation by removing the coolant 
temperature from the set of state variables. In reactivity feedback also, we 
consider the effect of only the reactor fuel temperature, and take the following 
equation intea d of Eq. (2.5) : 

p = - T° f ) (2.24) 

so that Eq. (2.10) is replaced by 

Xj = -bXj+bX 2 + aX 3 + aXjX 3 


(2.25) 



2.4.1. Model with Delayed Neutrons : 


The equations for this model are : 


X 1 = 

-bX j + bX 2 + aX^ + aX j X 3 

(2.26) 

X 2 = 

X 1 ' X 2 

(2.27) 

X 3 = 

cXj - cX 3 

(2.28) 

This is a three-dimensional continuous dynamical system with three 

parameters and two fixed points (0,0,0) and (-1, -1, -1). 

2.4.2. Model without Delayed Neutrons : 


As in Section 

Eq. (2.27), thus obtaining 

2.3.2, here also we put b = 0 in Eq. 

(2.26) and delete 

X 1 = 

aX 2 + aXjX 2 

(2.29) 

x 2 = 

cXj - cX 2 

(2.30) 


where we have written X 2 in place of of the previous system. It is under- 
stood that in the calculation of the parameter a, the value of 1 will be changed 
so as to partially compensate for ignoring delayed neutrons. 

This is a two-dimensional continuous dynamical system with two 
parameters and two fixed points (0,0) and (-1, -1). 

2.4.3. Model without Delayed Neutrons and with Constant Power Removal 
(Fast Transients) : 

In this model, it is assumed that heat from the reactor fuel elements 



is removed at a constant rate, equal to the equilibrium power. This results 


in dropping out of the X 2 term in Eq. (2.30), and we get the system : 


X 1 

= 

aX 2 + aXjX 2 

(2.29) 

*2 

= 

cXj 

(2.31) 


This system is two-dimensional and has only one fixed point i.e. 
(0,0), and two parameters. 

2.4.4. Higher-Order Feedback Models : 

For one-temperature higher -order feedback, we use the following 
equation for reactivity feedback (instead of the linear Eq. (2.24)) : 

P = a f (T f - T° ) + f(T f - T° ) (2.32) 

where f(x) is a nonlinear function with 

f(0) = f'(0) = 0. 

Now, instead of Eq. (2.25), we get 

Xj = -bXj + bX 2 + aX 3 + aXjX 3 + (1 + Xj) ^(X^ (2.33) 

where <j> is related to the function f in Eq. (2.32) and has the properties : 
4 >( 0 ) = 4 >'( 0 ) = 0 . 

The models with higher- order feedback have the same sets of 
parameters and fixed points as those with linear feedback. 
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2.5. Direct Feedback Models 

This model is a further simplification of the one-temperature feedback 
model. Here we remove the reactor temperature from the set of state variables 
so that a change in reactor power directly causes a feedback effect. To do 
this, we take the limit as 0 (very slow transients) in Eq. (2.23) and obtain 

Xj = X^. Substituting this in Eq. (2.25), we get 

X, = (a-b) Xj + bX 2 + aX 2 (2.34) 

2.5.1. Model with Delayed Neutrons : 

This model is given by 

Xj = (a-b) Xj + bX 2 + aX j (2.35) 

X 2 = Xj - X 2 (2.36) 

This system is two-dimensional with two parameters and two fixed 
points (0,0) and (-1,-1). 

2.5.2. Model without Delayed Neutrons s 

Neglecting delayed neutrons, we put b = 0 in Eq. (2.35) and ignore 
Eq. (2.36), thus obtaining 

X = aX + aX 2 (2.37) 

where we have written X in place of Xj. This is a one-dimensional continuous 
system with one parameter and two fixed points X = 0 and X = -I. 



2.5.3 Discrete Models : 


The two direct feedback models presented in this section will 
also be studied in their discretized versions. We can discretize the continuous 
systems by taking differences in place of differentials. Discretizing Eq. (2.37), 
we get 

x n+ l = (a 6 +1) X n + a 6 X 2 (2.38) 

where, 

6 = A x = step size of discretization, and 

X n » X n+ j = n th and n+1 th iterates of the variable X, respectively. 
This is a one-dimensional discrete dynamical system or map with one parameter. 

Similarly, we can discretize the system of Eqs. (2.35) and (2.36) 

to obtain 


X n+1 = 

(a 6 ■ 

- b 6 

+ 1) X + b 6 Y + a 6 X 2 
n n n 

(2.39) 

Y n+1 = 

6X n 

+ (1 

-6) Y„ 

(2.40) 


where, as before, 6 = A t is the step size. We have written X and Y for 
the variables Xj and X 2 of the continuous system and subscripts denote the 
number of iterations. This is a two-dimentional discrete dynamical system 
or map with three parameters. 



2.6. Values of the Parameters 


The dynamical systems, constructed in the previous sections involve 
some dimensionless parameters. These parameters have been defined in terms 
of the original physical parameters. The typical physical parameters for thermal 
and fast reactors are listed in Table 1. From the physical parameters, the 
dimensionless parameters may be calculated. 





CHAPTER 3 


INVESTIGATIONS ON THE DYNAMICAL SYSTEMS 
FOR NUCLEAR FISSION REACTORS 

In the previous chapter, we constructed dynamical systems ranging 
from one-dimensional to four-dimensional continuous and discrete systems. 
Now we proceed to investigate the behavior of these systems by analytical 
as well as numerical methods, with an eye on the possibility of the existence 
of limit cycles. The requisite concepts for analysis and computation have been 
taken from Refs. [14-27] and briefly mentioned in Appendix-A. 

3.1. One-Temperature Feedback Models 

First of all, we study the dynamical systems of Section 2.4. 

3.1.1. Model without Delayed Neutrons : 


This is a two-dimensional system of equations 


*1 

= 

aX 2 + aXjX 2 

(3.1) 

*2 

= 

cXj - cX 2 

(3.2) 


having two fixed points (0,0) and (-1, -1). 
Local Behavior at the Operating Point (0,0) 


The Jacobian matrix at this point is 
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whose characteristic equation is 

X 2 + cX - ac = 0 (3-3) 

having the roots (eigenvalues) 

X = - c/2 ± ( 1 / 2 ) IT 2 ~~~ 

l >2 'J c + 4 ac 

From this, we can infer that 

(1) This point is unstable for a > 0 (assuming c > 0) and locally stable 
for a < 0. 

(2) From the movement of the eigenvalues in the complex -plane [Fig.l], 
we see that for a > -c/4 (and c > 0), Xj and X 2 are real. At a = -c/4, 
they merge and then move vertically away from the real axis 
as 'a' decreases further. Therefore, locally (i.e. for small disturbances) 
the operating point behaves like an unstable node for a >0, a stable 
node for -c/4 < a <0, and a stable focus for a < -c/4 with oscillations 
becoming faster as'a' decreases, 

2 

(3) To get pure imaginary eigenvalues, we must have c < 4ac and 
c = 0, which is not possible. Therefore, Hopf bifurcation as well 
as existence of any periodic orbits in the neighborhood of the origin 
(i.e. the operating point) is ruled out. 

Local behavior at the Shutdown Point (-1, -1) 

The Jacobian matrix at this point has the eigenvalues 

Aj = - a and X 2 = - c. 



Therefore, 


(]) This point is locally stable for a > 0 and unstable for a < 0. 

(2) Since the eigenvalues are real, no Hopf bifurcation is possible 
and periodic orbits in the neighborhood of the shutdown point are 
rule out. 

Global Behavior 

If this system has no limit cycles, the two fixed points constitute 
the global set of attractors and repellers. By combining the local behavior 
at the two fixed points, we can draw the global phase-portrait. 

Numerical Experiments 

To verify this global phase-portrait, we performed numerical experi- 
ments with a = -1.0 and c = 1.6 taking several initial conditions. The resulting 
plot is shown in Fig. 2. It confirms the existence of a stable focus at the 
operating point (0,0,0). For simulating the dynamical system, IM5L routine 
DVERK based on Verner's method was used with an interactive main program. 

Some Further Remarks 

Although we found no Hopf bifurcation in the system of Eqs. (3.1) 
to (3.2), we can create a Hopf bifurcation by slightly modifying the system 



where d is a new parameter, d = 0 corresponds to the original system. The 
characteristic equation of the new system at (0,0) is 

(d-c) ^ - ac = 0 (3.5) 


The eigenvalues now are 



= << 1 - 0/2 * 0 / 2 ) / (d . c) 2 + , ac 


which become purely imaginary at d = c, causing a Hopf bifurcation (assuming 
a<0 and c > 0, which represents a realistic case). As d increases beyond c, 
the eigenvalues go from left-half to right-half of the complex-plane as shown 
in Fig. 3. Therefore, from the Hopf bifurcation theorem, it follows that there 
exists either a stable limit cycle for d > c or unstable limit cycle for d < c 
or a nonlinear center at d = c. However as we will see in the next section, 
we find periodic solutions for the case d = c. It is on this basis that periodic 
solutions in the system of Eqs. (3.1) - (3.2) are ruled out globally. 


3.1.2. Model without Delayed Neutrons and with Constant Power Removal 
This model is given by the system of equations 

Xj = aX 2 +aXjX 2 (3.6) 

X 2 . = cX, (3.7) 

which has only one fixed point (0,0). We note that this system can also be 
obtained from the system of Eqs. (3.1) and (3.4) by setting d = c for which 
we indicated the possibility of a nonlinear center. At (0,0), this system has 



the Jacobian matrix 

' 0 a 
c 0 

whose characteristic equation is 

X 2 - ac = 0 (3.8) 

and the eigenvalues are 

X j ^ — i / a c 

which are purely imaginary for a < 0, c > 0. Therefore, the linear system has 
a center at (0,0) which may or may not be retained in the nonlinear system. 

Numerical Experiments 

To find out whether or not a nonlinear center exists, we performed 
numerical experiments with a = -1, c = 1.6 taking several initial conditions 
and each time found a closed curve around (0,0), thus indicating the existence 
of a nonlinear center [Fig. 4]. For some initial conditions, the curve went below 
X 2 = -1 line which is unrealistic. 

3.1.3. Model with Delayed Neutrons : 


Now 

we come to a three-dimensional system 

of the equations 

s 

X 1 

= -bX. + bX_ + aX, + aX,X 0 

1 z i J 5 

(3.9) 

X 2 

(N 

X 

1 

xT 

ii 

(3.10) 

•X 

cX ] - cX 3 

(3.11) 


having two fixed points (0,0,0) and (-1, -1, -1). 
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Local Behavior at the Operating Point (0,0,0) 

* 

The Jacobian matrix at this point is 

-b b a 
1 -1 0 
_ c 0 -c 

whose characteristic equation is 

X^ + (c + b+l)A^ + c(b+l-a)A-ac = o (3. 12) 

For studying the movement of eigenvalues in the complex -plane, we rearrange 
the above equation as 

A (A + c) ( A + b + 1) = ac (A + 1) (3.13) 

For the case a = 0 (no feedback), we find that all the three eigen- 
values are real 

A j = 0, A 2 = - c, A 3 = - b - 1 

For the other cases, a < 0 (negative feedback) and a >0 (positive 
feedback) we can find the prohibited regions for the eigenvalues on the real 
line and hence investigate what happens as a decreases or increases starting 
from zero. We identify four cases : 

(1) Case a<0, c > 1 (assuming b + 1 > c) : 

A <£(-<= 0 , -b -1] U[ -c, -1] U[0, co ) (3.14) 

From this restriction and with the help of numerical experiments, 
we find that, as a decreases from a = 0, Aj tends to -1 ; whereas A 2 and A^ 
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approach each other, merge at a negative value and then move vertically 
away from the real axis [Fig. 5a]. Therefore, the operating point behaves locally 
as a stable node or a stable focus. 

(2) Case a < 0, 0 < c < 1 (assuming b > 0) : 

X <£ (- 00 , _ b - 1] U [-1, -c] U [0,°° ) (3.15) 

From this constraint and with the help of numerical experiments, 
we find an interesting movement of eigenvalues in the complex-plane [Fig.5.b]. 
As a decreases from a = 0, A j and A^ approach each other, merge, become 
complex conjugates, move to the left and land again on the real axis in the 
interval (-b-1, -1). Then one of the eigenvalues tends to -1, whereas the other 
merges with the third eigenvalue and forms complex conjugates which move 
away from the real axis vertically. Therefore, the operating point locally behaves 
as a stable node or a stable focus. 

(3) Case a > 0, c >1 (assuming b+1 > c) : 

X f[-b-l, -c] U [-1, 0] (3.16) 

From this prohibition we may conclue that no two eigenvalues 
are free to merge and form complex conjugates, and one of the eigenvalues 
is positive. Therefore, the operating point locally behaves like an unstable 
node. 

(4) Case a > 0, 0 <c < 1 (assuming b > 0) : 


X ^ [-b-1 , -I] U[-c, 0] 


(3.17) 



From this restriction, we may conclude, as in the previous case, 
that no two eigenvalues are free to merge and form complex conjugates and 
one of the eigenvalues is positive. Therefore, the operating point locally behaves 
like an unstable node. 

Local Behavior at the Shutdown Point (-1, -1, -I) 

The Jacobian matrix at this point is 

-b-a b 0 
1 -1 0 
c 0 - c 


whose characteristic equation is 



( A + 

c) [A + (a+b+l)A + a ] = 0 

(3.18) 

having the eigenvalues 



S 

= 

- c, 


*2,3 

= 

- (a+b+l)/2 ± (1/2) y (a+ b +l) 2 _ ttQ 



= 

- <a + b*l)/2 ± (1/2) / (a+b _, )2 + „ b 



From physical considerations, we must have b > 0. Therefore, all 
the eigenvalues are real and all of them are negative for a > 0, but not for a < 0. 
Therefore, locally, the shutdown point behaves like a stable node for a >0 
and an unstable node for a < 0. 

Hopf Bifurcation at the Shutdown Point 

Pure imaginary eigenvalues at (-1, -1, -1) occur at b = -a-1 for 
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a > 0 and the derivative 

Re ( i2 ) < 0 

Therefore, for fixed a >0, c >0, this is a case of Hopf bifurcation in which 
a pair of complex conjugate eigenvalues crosses the imaginary exis from right 
to left as b increases beyond the critical value b + = -a-1 while the third eigen- 
value remains in the left-half plane . From Hopf bifurcation theorem, 

it follows that there exists either a stable limit cycle for b >b + , or an unstable 
limit cycle for b < b + or a nonlinear center at b = b + . We are interested 
only in b > 0 for which a stable limit cycle may or may not exist. Even if 
it exists, it is likely to be unrealistic, because it would surround the shutdown 
point (-1, -1, -1), sometimes going below it. 

Hopf Bifurcation at the Operating Point 

To locate Hopf bifurcation at (0,0,0), we suppose Eq. (3.12) has 
the roots A . _ = a ± i 3 » A , = r with r, a , 3 real. These roots are related 

1 ,2 j 

by 

Aj + A 2 + ^3 = ~ (c+.b+l) 

Aj * 2 + X 2 X 3 + X 3 X I = C (b+1 ' a) 

X 1 X 2 X 3 = aC 

That is, 

2a + r = - (c+b+1) (3.19) 

= c (b+l-a) 


a 2 + 3 2 + 2ar 


( 3 . 20 ) 



2 

For pure imaginary eigenvalues, a + = 0 and 8 + > 0, therefore 


r + • = -(c+b+1) 

(3.22) 

2 

g + = c (b+ 1 -a + ) 

(3.23) 

2 

6 r = a c 

(3.24) 

Therefore, 


c(b+l-a + ) (c+b+1) = a + c 


or 


a + = (b+1) (c+b+1 )/(c+b) 

(3.25) 

2 

This is the Hopf bifurcation point provided that g > 0, i.e. 



c(b+ l-a + ) >0 
or 

ca + < c(b+l) 
or 

c(b+l) (c+b+1 )/(c+b) < c(b+l) 
Assuming that c > 0 and b > 0, we must have 

c+b+1 < c+b 
i.e. 1 < 0 


which is impossible. 



Therefore Hopf bifurcation as well as the existence of any periodic 
orbits in the neighborhood of the operating point is ruled out. 

Global Behavior 

For negative feedback (a < 0), no limit cycles are expected surround- 
ing either of the fixed points. Therefore, we can combine the local behavior 
at the two fixed points to get the global phase-portrait. It will be a combination 
of a stable node or a stable focus at (0,0,0) and an unstable node at (-1 ,— 1 , 1 ). 

For positive feedback (a > 0), there is no limit cycle around (0,0,0), 
but there may be a stable limit cycle around (-1, -1, -1). Therefore, disturbances 
from the operating point may either become unbounded or asymptotically 
tend to the limit cycle. 

Numerical Experiments 

Numerical simulations were performed to verify the phase-portraits 
predicted above. The parameter values were taken for a thermal reactor 
(Chooz Sena PWR) : a = -2614, b = 2143, c = 1.678 taking several initial 
conditions. The resulting orbits are given in Table 2 and their projection on 
Xj - X ^ plane is shown in Fig.6. The operating point behaves like a globally 
stable node, attracting all orbits. 

3.1.4. Higher-Order Feedback Models : 

The dynamical systems obtained from these models are the same 
as those for the corresponding linear feedback models but for an additional 



term in the first equation. The inclusion of this additional nonlinear term, 
from which the linear part was already separated, does not affect the eigenvalues 
at the operating point (0,0,0) but may affect those at the shutdown point 
(~1> > ~J)* Therefore, the local behavior at (0,0,0) is not affected but that 

at (-1, -1, -1) may be affected. The Hopf bifurcation at (0,0,0) will occur 
at the same parameter values as for the linear feedback models but at (-1, -1, -1), 
the Hopf bifurcation point may change. Moreover, the direction of branching 
and the stability of limit cycles may change due to the inclusion of higher 
order feedback. These things may affect the global behavior. 

3.2. Two-Temperature Feedback Models 

Now we study the dynamical systems of Section 2.3. 

3.2.1. Model without Delay Neutrons : 


This is a three-dimensional system of equations 


X 1 

= aX^ + qX^ + aXjX^ + qXjX^ 

(3.26) 

*2 

= (k-h) Xj - kX 2 + hX 3 

(3.27) 

X 3 

= fx 2 - fx 3 

(3.28) 

having two fixed points (0,0,0) and (-1, -1, -!)• 


Local Behavior at 

the Operating Point (0,0,0) 



The Jacobian matrix at this point is 



whose characteristic equation is 


X 3 + (k+f) X 2 + (f-a) (k-h) X - ( a +q) (k-h) f = 0 


(3.29) 


For the case a + q = 0, we get the eigenvalues 


X 


1 




- - (k+f)/ 2 ± (1/2) 


+ 4f + 4a (k-h) 


Since f > 0, h > 0 and k-h > 0, we have real eigenvalues for 


a > - 


(f-k) 2 + 4hf 
4 (k-h) 


= -a (say) 


For a < - a , we have a complex conjugate pair of eigenvalues with negative 
real part, which move vertically toward the real axis as a increases, merge 
and then move along the real axis in opposite directions. One of them becomes 
positive for a > f . All the time, the eigenvalue ^ remains zero [Fig. 7], 

Thus, we can say that for a+q = 0, the local behavior of the 
system around (0,0,0) cannot be determined by linear analysis only (since one 
of the eigenvalues is zero). 


Hopf Bifurcation at Operating Point 


Since k>0 and f > 0, no pure imaginary eigenvalues are possible 
for the case a+q = 0. Therefore, there is no Hopf bifurcation in this case. 

Suppose we fix a+q = G = a constant and define new parameters 
(constants) F = k+f >0 and H = k-h > 0. Then Eq. (3.29) becomes 



whose roots A,, A 2 , A3 are related by 


Aj + A 2 + A 3 = - F 


Aj A 2 + X 2 A3 + A 3 Aj = H (f-a) 
Aj A 2 A 3 = GHf 


(3.31) 


(3.32) 


(3.33) 


In the vicinity of Hopf bifurcation point, we have a pair of complex conjugate 

eigenvalues, say Aj ^ 2 = a± ig , g 2 > 0 and a real eigenvalue, say a 3 = r. Then 
we have 


2 a + r = - F 


a 2 + g 2 + 2 a r = H (f-a) 


( a 2 + g 2 ) r = GHf 


(3.34) 

(3.35) 

(3.36) 


At a = a + (Hopf bifurcation point), we have a + = 0, g + = to + > 0. Putting 
these in the above equations and eliminating r + and cd + , we get 


cu 


f (G+F) 


GHf 


Since H, F and f are positive, we must keep G < 0 to have w + > 0. We also 
note that the real eigenvalue is negative at the bifurcation point. 

To determine the direction of movement of the eigenvalues in the 
complex-plane, we differentiate Eqs. (3.34) - (3.36) with respect to a, and 
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N3 

+ 

II 

0 

(3.37) 

2a a' + 2 6 6'+ 2 a r’ + 2r a ' = - H 

(3.38) 

r (2 aa' + 2 66') + ( a 2 + $ 2 ) r' =0 

(3.39) 


Putting the values of a, 3 > r at a = a + , and solving for a’ we get 

a' (a + ) = -^( 1 - ) _I (3.40) 

which is positive for H, F, f >0 and G < 0. Therefore, the complex conjugate 
eigenvalues cross the imaginary axis from left to right and the real eigenvalue 
remains in the left-half plane as a increases beyond a + . From Hopf bifurcation 
theorem it follows that there exists either a stable limit cycle for a > a + 
or an unstable limit cycle for a <a + or a nonlinear center at a = a + . [Fig. 8]. 

Numerical Experiments 

Numerical simulations were performed for G = -5, k = 2, h = 0.5, 
and f = 8 for which a + = 4. Trajectories were computed taking a = 4.5 and 
two initial conditions. The resulting orbits are given in Table 3,4 and projection 
on Xj - X 2 plane is shown in Figs. 9,10. The orbit starting outside the limit 
cycle spiralled into it, and the one starting inside the limit cycle also spiralled 
out to it. The dimensionless times taken by these orbits to practically approach 
the limit cycle were 80 and 120 respectively. The period of the limit cycle 
was approximately 12. This indicates the existence of a stable limit cycle 
for a > a + . 



39 


3.2.2. Model with Delayed Neutrons j 

This is a four-dimensional dynamical system of the equations 
Xj = -bXj + bX 2 + aX 3 + qX^ + aXjX 3 + qXjX^ 
X- = X, - X, 


(k-h) Xj - kX 3 + hX^ 


fx 3 - fX^ 


having two fixed points (0,0, 0,0) and (-1, -1, -1, -1). 


Characteristic Polynomial at the Operating Point (0,0,0, 0) 


( 3 . 41 ) 

(3.42) 

(3.43) 

(3.44) 


The Jacobian matrix at this point is 


-b b a q 

1-10 0 

k-h 0 -k h 

0 0 f -f 


whose characteristic polynomial is given by 


where 


P(X) = A 4 + a'j X + a 2 X + a 3 X + a^ 


1 + b + k + f 
f + k (1+b) + (a+f) (k-h) 
a + f(l+b+a+q) 

(a+q) (k-h) f 


(3.45) 



Locating Hopf Bifurcation by Direct Decomposition Technique 

In this technique, we decompose the characteristic polynomial near 
the Hopf bifurcation point as 

P(A) = ( A 2 + u) ( A 2 + p A + p 2 ) + A A + B (3.46) 

coefficients in this equation are given recursively by 


P-i = °» p c = 1 ; 

Pj^ = a^ - w Pj < _2 ’ k = 1 , 2, .... , n-2 ; 

A = a„ , - wp n , ; B = a - um - 
n-1 r n-3 n mi-2 


In our case, n = 4, and 


P-1 = 0 s Po = 1 

Pj = aj - up j = 3j (3.47) 

p 2 = a 2 - a) p o = a 2 - a) (3.48) 

A =3^-( J )P^=a^-u)a i ( 3.49) 

B = - up 2 = a^-a)(a 2 -w) 

= a^ - a 2 aj + to (3.50) 


At the Hopf bifurcation point, the parameters should be such that A = 0, 
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From Eq. (3.51) we get 

w = a^/a j 

A 

Substituting the value of w in Eq. (3.52) from Eq. (3.53), we get 
a^ - (a^/aj) (a 2 - ayaj) = 0 


or, 


2 2 
a l a 4 ' a l a 2 3 3 + a 3 


0 


or, 


a l a 2 a 3 “ 


a l a 4 


+ a 


2 

3 


Moreover, for pure imaginary roots, we must have u > 0 
i.e. , 

a^aj > 0 


(3.53) 


(3.54) 


(3.55) 


Eqs. (3.54) and (3.55) are the necessary conditions for Hopf bifurcation in the 
system of Eqs. (3.41) - (3.44). 

Thus, we see that there is a possibility of Hopf bifurcation at a = a . 
Periodic orbits around the operating point, therefore, cannot be ruled out in 
the vicinity of the operating point. Further analysis, however, appears to be 

too involved. 

3.2.3. Higher-Order Feedback Models : 

In the two-temperature feedback also, nonlinear feedback models 
may be considered. Here too, the effect of inclusion of higher-order feedback 



will be similar to that in the one-temperature feedback case. The location 
of Hopf bifurcation point will remain unchanged but the direction of branching 
and stability of the limit cycles may change. Further bifurcations from limit 
cycle may also be affected by the extra nonlinear terms. 

3.3. Direct Feedback Models 

Now we study the simplest models i.e. those of Section 2.4. 

3.3.1. Model without Delayed Neutrons : 

This is the simplest of the models, as given by the equation : 

X = aX + aX 2 (3.56) 

having two fixed points X = 0 and X = -1. Its Jacobian is 

dX/dX= a at X = 0 

dX/dX= -a at X = -1 

Therefore, X = 0 is unstable for a > 0 and locally stable for a < 0. Except for 
unrealistic initial conditions, X < -1, it is also globally asymptotically stable 
for a < 0. At a = 0, an exchange of stability occurs between the two fixed 
points. 

Since this systems is a one-dimensional continuous dynamical system, 
it cannot have any limit cycles or periodic orbits. 
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3.3.2. Model with Delayed Neutrons : 

This is a two-dimensional system of equations 

= (a-b) Xj + bX 2 + aX* (3.57) 

X 2 = Xj - X 2 • (3.58 ) 

having two fixed points (0,0) and (-], -1) . 

At the operating point (0,0) the Jacobian matrix is 

[•? -i] 

whose characteristic equation is 

X 2 + (b-a-1) \ - a = 0 (3.59) 

having the eigenvalues 

*1,2 = -<'« (b-a-D * (1 «/ (bs5+1) 2 . « a 

= -0/2) (b-a-1) - <l/2)/ (b-a-l)^ 7 4b 

For b > 0, both the eigenvalues are real. In addition, if a <0, they are 
negative, hence a stable node. Keeping b >0, as a increases past a = 0, one 
eigenvalue becomes positive. Therefore, there is a saddle-node bifurcation at 


a = 0. 



Hopf Bifurcation at the Operating Point 
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Pure imaginary eigenvalues occur at b = a-1 for a fixed a < 0. As 
b increases beyond this value, the eigenvalues cross the imaginary axis from 
right to left. Therefore, from Hopf bifurcation theorem there exists either 
a stable limit cycle for b <0 or an unstable limit cycle for b > 0 or a non- 
linear center at b = 0. We are interested in b > 0 for which an unstable limit 
cycle is possible. 

3.3.3. Discrete Models 


The discrete model without delayed neutrons is the one-dimensional 


map 


X n+] = (a 6 + i) X n + a 6 X x 


which can be transformed to the well known logistic map 


x n+l = 11 x n " x n> 


by the transformation 


* < •air i 1 > x 


where 


a 6 + 1 . 


(3.60) 


(3.61) 



The logistic map is known to have chaotic behavior for y > 4. This corresponds 
to a 6 >3, which means that a must be positive, and the dimensionless time 
step, 6 , should be of the order 1 (because a is of the order 1), which is unrealis- 
tic discretization. 
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The discrete model with delayed neutrons is the two-dimensional 

map 

X n+1 = (a 6 - b$ +l)X n +b6 Y n + a5 X* <3.«> 

Y n + 1 = 6X n + 6 > Y n (3.63) 

which can be transformed to the Henon map (also chaotic), taking 6 = I, 

x n+ 1 = l + y n + oxJ (3.64) 

I'n.l = bx n < 3 - 6 » 

by the transformation 

x = AX - B, 

y = b (AY - B) 

where 

A = a/a , B = 2/(b-l), 

■ i ■ 

a = (1/4) (b- 1 ) (b-a- 1 ) 

Here we had to take the dimensionless time step, 5 , equal to unity, which 

is an unrealistic discretization. 

Thus, we see that in both the discrete models discussed above, 
we had to take unrealistic discretizations to see chaos. The discrete models 
based on discrete population dynamics may exhibit different kind of behavior. 


* * * * 



CHAPTER 4 


SOME ADDITIONAL INVESTIGATIONS 

In addition to analysing the fission reactor dynamical systems 
for periodic behavior, we have, during the course of this work, also attempted 
to see if these systems, have any chaotic behavior. Methods for identifying 
chaos were studied. The method of Lyapunov exponents was applied to the 
two-temperature feedback model without delayed neutrons. 

4.1. Methods for Identifying Chaos 

In the last decade, new phenomena have been observed in all areas 
of nonlinear dynamics, principal among these being chaotic oscillations. Chao- 
tic oscillations are the emergence of bounded, aperiodic random like motions 
from completely deterministic systems. Such motions had been known in fluid 
mechanics, but they have recently been observed in low-order mechanical 
and electrical systems and even in simple one-degree-of-freedom problems. 
The random character of the motions in these deterministic systems arises 
from sensitive dependence on initial conditions. This leads to exponential diver- 
gence of nearby trajectories. 

Methods for identifying chaos include Lyapunov exponents, Poincare 
map, Fourier spectrum, and the identification of routes to chaos through bifur- 
cations. Details of these methods may be found in Refs.[21-24, 28-30]. 
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4.1.1. Lyapunov Exponents : 

Numbers that measure the exponential attraction or separation 
in time of two adjacent trajectories in phase space with different initial condi- 
tions are called Lyapunov exponents. 

If one imagines a set of initial conditions within a sphere of radius e in 
phase space, then for chaotic motions, trajectories originating in the sphere 
will map the sphere into an ellipsoid whose major axis grows exponentially 
with an exponent X > 0, which is known as a Lyapunov exponent. If d Q is a 
measure of distance between two starting points, at a small but later time 

d (t) = d o 2 Xt (4.1) 

The divergence of chaotic orbits can only be locally exponential, 
since if the system is bounded, as most physical experiments are, d (t) cannot 
go to infinity. Thus, to define a measure of this divergence of orbits, we must 
average the exponential growth at many points along a trajectory. 

4.1.2. Poincare' Sections : 

A Poincare' section is the sequence of points in phase-space generated 
by the penetration of a continuous evolution trajectory through a generalized 
surface or plane in the space. For a 3-dimensional case, we can take a plane, S, 
and place the points of intersection corresponding to a given direction of 
crossing. The plane is chosen so that the trajectory T intersects S at P Q , P,,.~ , 



where the dynamics of the system are assumed to be such that the trajectory 
continually crosses from one side of S to the other. S can be any plane, but 
an appropriate choice yields sections that are more easily analysed. In a three- 
dimensional phase space, the Poincare 1 section of a periodic orbit is a single 
point or a set of a finite number of points. For a queasiperiodic orbit, it is a 
continuous closed curve, but is not traversed continuously by successive inter- 
section points of the trajectory with the plane. For an aperiodic orbit, the 
Poincare' section is a set of points distributed along an unclosed curve. It 
can be reduced further to a one-dimensional first return map by defining a 
coordinate for each point on the curve and plotting each point of intersection 
against the previous one. 


fr.1.3. Fourier Transform and Power Sepctrum : 


The Fourier transform of a discrete time series x. is defined to 
be the operation creating a corresponding discrete series such that 

x^ = ZZT j^] x j ex P ( ^ ) > k = 1, ... , n. (4.2) 

The power spectrum is the graph of as a function of frequency f (f=k. A f) 
where is defined by 


2tr mk 


Z i / ZH 1 1 \ 

z = l . ^ COS ( ) 

m= m n 


(4.3) 


where is the autocorrelation function of the signal x. defined by 
r m 1 



For periodic signals (with period T) we get a peak at the frequency 
l/T, its sidelobes, and possibly a certain number of other peaks (and their 
sidelobes) that are harmonics of the fundamental frequency. Fora quasiperiodic 
signal (with fundamental frequencies f } and y several peaks at independent 
frequencies, f = m^ + rr^y where nrij and m ? are positive integers. If f J /f 2 
is irrational, then the power spectrum is dense. For an aperiodic signal, the 
power spectrum is continuous. But this is not a conclusive evidence of an 
aperiodic or nonperiodic signal, because the power spectrum of a quasiperiodic 
signal with very high number of frequencies has a similar appearance. 

fr-2. Calculation of the Largest Lyapunov Exponent 

The signs of the Lyapunov exponents provide a qualitative picture 
of a system's dynamics. In a three-dimensional continuous dissipative system 
the only possible spectra, and the attractors they describe, are as follows 
: (+, 0, -), a strnage attractor (chaos) ; (0, 0, -), a two-torus (quasiperiodic); 
(0, -, -), a limit cycle; and (-, t, -), a fixed point. Since the largest Lyapunov 
exponent dominates the behavior of a dynamical system, its calculation is 
important to determine the nature of the system. If it turns out to be positive, 
it confirms that the system exhibits chaos. 

Suppose we want to calculate the largest Lyapunov exponent of 

a flow s 

x = F (x) 


whose variational equations are : 



A 6 x 


(4.6) 


6x = 


where A is the Jacobian matrix. To calculate the largest positive Lyapunov 
exponent, we have to measure the divergence of nearby trajectories averaged 
over a long time. Eq. (4.6) is to be integrated along a trajectory over a small 
enough time interval starting with arbitrary initial conditions. Redefining the 
initial conditions, the process is to be continued along the same trajectory. 
Taking a sufficiently large number, N, of such steps, the largest Lyapunov 
exponent may be calculated as 


X 



v d( V 

k=l log 2 d(tj^j) 


(4.7) 


The methods for calculating the complete spectrum of Lyapunov 
exponents of a system are discussed in Refs. [28-30]. 


4.2.1. Algorithm and Program : 

Based on the method give above, an algorithm for calculating the 
largest Lyapunov exponent was written which is as follows : 

1. Starting from a. point in the state-space, calculate the trajectory 
over a sufficient time to allow the orbit to approach the attractor. 

2. Choose an arbitrary variation vector of unit magnitude. 

3. Starting from the present point in the state-space, integrate the 
equations of the flow and its variation simultaneously, over a suitably 


chosen time step. 



ACC. No. r’k, 

4. Calculate the magnitude, d (t,), of the variation vector, where 

K 

k is the present step number. 

5. Calculate log 2 d (t^) and add it to the sum of all the previous 
values of d(t), i = 1, ..., k. 

6. Divide the new sum by the time over which the integration has 
been done to obtain the latest estimate of the largest Lyapunov expo- 
nent. 

7. Normalize the variation vector to unit magnitude. 

8. Repeat steps 3 to 7 until a satisfactory estimate of largest Lyapunov 
exponent is obtained. 
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An interactive program [Appendix-B] based on the above algorithm 
was written in FORTRAN IV and implemented on DEC 1090 computer. The 
program was tested for the Lorenz and Rossler systems (both chaotic), and 
the results agreed with those in Refs. [28, 29]. 

The calculation of the complete spectrum of the Lyapunove exponents 
is described in Refs. [28-30]. The program of Wolf et al [28] was made inter- 
active and implemented on DEC 1090 computer for the Lorenz and Rossler 
systems. The results agreed with those of the author. 

4 . 2 . 2 . Calculation of the largest Lyapunov Exponent for the Two-Temperature 
Feedback Model without Delayed Neutrons . 

For this reactor model, Hopf bifurcation had been shown analytically 
in Section 3.2.1. and the existence of a stable limit cycle observed by numerical 



Lyapunov exponent for this system taking two sets of values of parameters. 
Tim parameter values were chosen such that the operating point would become 

unstable and Hopf bifurcation theorem would predict the existence of a stable 
limit cycle. 

I . Parameter Set No. 1 

G = -2, k = 2 , h = 1.6, f = 10.0 for which Hopf bifurcation occurs 
ata + =8-|j ; a = 9 was chosen. 

First Attempt 

-3 

Execution was done for 80,000 time steps of size 10 
skipping the first 1000 steps and printing every 100th step. The tolerance 
chosen to be 10 . The execution on DEC 1090 took 5.24s of CPU time. 

The results of this attempt are shown in Table . Towards 
the end of the execution, the estimate of the Lyapunov exponent was seen 
to be oscillating between positive values only (+0.02 and +0.10). It was suspected 
that, for these parameter values, further bifurcations of the limit cycle had 
taken place, and the existence of a chaotic attractor was suspected. To confirm 
this, another attempt with a longer time of evolution was made. 

Second Attempt 

This time, 6000 steps of size 1 were taken, skipping the 
first 100 steps and printing every 100th step. The tolerance was chosen as 
I0 - ^. The execution took 2 min. 11.24s of CPU time on DEC 1090. 
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Towards the end of the execution, the estimate of the largest 
Lyapunov exponent was seen to oe oscillating between small positive and nega- 
tive values, of the order 10 , which can be taken to be approaching zero 

Hence, the existence of a chaotic orbit is ruled out and that of a limit cycle 
or a torus is expected. This demonstrates that, taking insufficient time of 

evolution may lead to erroneous conclusions. 

2. Parameter Set No.2 

G = -5, h = 0.5, k = 2, f = 8 for which Hopf bifurcation occurs 
at a + = 4 ; a = 4.5 was chosen. 

\ 

For this parameter set, numerical experiments had already been 
performed and reported in Section 3.2.1, where a stable limit cycle had been 
observed. The execution for calculating the largest Lyapunov exponent was 
done for 4000 time steps of size 1, skipping first 100 steps and printing every 

100th step. As for the previous parameter set, here also the estimate was 

seen to be oscillating between small positive and negative values of the order 
I0"\ which can be taken to be approaching zero. Hence, the existence of 
a chaotic orbit is ruled out in this case also, and a limit cycle or a torus 
is expected. This, in conjunction with the results of numerical simulations, 
may be used to infer that a stable limit cycle exists for this set of parameter 


values. 



CHAPTER 5 


CONCLUSION 

In the preceeding chapters nuclear fission reactors were modelled 
using point reactor kinetics and reactivity feedback effects. Transforming 
the equations to dimensionless form, dynamical systems for various reactor 
models were constructed. These dynamical systems were then investigated 
for the possibility of the existence of limit cycles by analytical methods as 
well as numerical experiments. In one of the models (the two-temperature 
feedback model without delayed neutrons) a stable limit cycle was predicted 
using the method of Hopf bifurcation. This was confirmed by numerical experi- 
ments. Some additional investigations on the possibility of chaos were also 
done, by computing the largest Lyapunov exponent. The result was in the negative 
for the parameter values chosen (after the Hopf bifurcation). 

5.1. Summary of the Dynamic Modeling 

The dynamic models were obtained using the equations of point 
reactor krnetics and energy balances. Tne feedback effects of the energy balance 

equations were included in the point kinetic equations. 

«c thnc obtained were transformed in such 
The dynamic equations thus oDtameu wv 

a way that all the state variables, the independent variable, and the parameters 

became drmens.onless, while still allowing direct physical interpretation. The 
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state variables were normalized with respect to the equilibrium operating 
point of the final system (after the autonomous disturbances) and this point 
was translated to the origin. The time was normalized with respect to the 
decay constant of delayed neutron precursors. Thus, the number of parameters 
in the system was greatly reduced. The shutdown point was also made indepen- 
dent of all the parameters by choosing the normalizations appropriately. 

5 . 2. Summary of the Investigations 

The dynamical systems constructed for various reactor models 
were investigated by analytical methods as well as numerical experiments. 
The movement of eigenvalues was studied with the aim of locating a Hopf 
bifurcation. The direct decomposition technique was used to locate Hopf bifur- 
cation in the two-temperature feedback model with delayed neutrons. Numerical 
experiments were performed to verify the theoretical predictions. The results 
are as follows : 

]. One-Temperature Feedback Models 

a) without delayed neutrons : Hopf bifurcation was ruled out. A stable 
focus at the origin was predicted and experimentally confirmed. 

b) without delayed neutrons and with constant power removal : Possibility 
of a nonlinear center was indicated and experimentally confirmed. 

c) with delayed neutrons : Hopf bifurcation was ruled out. A stable 
node at the origin was ' predicted and experimentally confirmed. 



d) higher order feedback models : Compared with the linear feedback 
models. 

2. Two-Temperature Feedback Models 

a) without delayed neutrons. Hopf bifurcation was located and periodic 
orbits were predicted. Numerical experiments indicated the occurrence 
of a supercritical Hopf bifurcation and the existence of a stable 
limit cycle. 

b) with delayed neutrons : The necessary condition for the Hopf bifur- 
cation was determined by using the direct decomposition technique. 

c) higher-order feedback models : Compared with the linear feedback 
models. 

3. Direct Feedback Models 

a) without delayed neutrons : One-dimensional system, hence no limit 
cycles are possible. 

b) with delayed neutrons : Hopf bifurcation was located and the existence 
of an unstable limit cycle was predicted. 

c) discrete models : the above two models were discretized and reduced 
to the logistic map and the Henon map respectively, both of which 
are well known for exhibiting chaos. However, it was noted that, 
in our systems, chaos can occur only for unrealistic discretization. 


f 


4. For detecting chaos, in the continuous systems, a program for computing 
the largest Lyapunov exponent was written and tested. The computations for 
the two-temperature feedback model without delayed neutrons yielded a zero 
value of the largest Lyapunov exponent, thus ruling out the possibility of chaos 
for the parameter values chosen, and further confirming the existence of 
a stable limit cycle. 

5.3. Suggestions 

For further investigations in the direction of the present work, 
the following suggestions can be given : 

1) The direction of branching, the stability of the limit cycles and 
their domains of attraction need to be investigated further. 

2) Absence of Hopf bifurcation does not exclude periodic orbits of 
large amplitude and periods. This needs further investigation. 

3) The possibility of chaos for realistic values of the parameters 
may be investigated by computing the largest Lyapunov exponent. 
If one had access to a supercomputer, • one could calculate the 
largest Lyapunov exponent for many values of the parameters 
and thus determine the ranges of parameter values for which chaos 
exists. 

4) The nature of the attractor may be studied by taking a Poincare 
section. 

Si The behavior of the dynamical systems may also be investigated 
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by taking the Fourier transform and studying the power spectrum. 

5Jf Epilogue 

In the present work, several mathematical models based on different 
physical assumptions have been constructed and studied, with the objective 
of locating limit cycles. Different models exhibited different kind of behavior, 
ranging from globally stable fixed points to centers and limit cycles. Further 
investigations on mathematical models may be carried out, but it is important 
to note that they are based on certain idealizations and approximations. Predic- 
tions based on mathematical models need to be verified by experiments on 
the actual physical systems. 


* * * * 



APPENDIX - A 


SOME MATHEMATICAL CONCEPTS 
A.I. The Hopi Bifurcation Theorem 

Let us consider a system of ordinary differential equations 

= F (x,y ) (A.I) 

which has a stationary point x = 0 € R n and the critical value of the bifurcation 
parameter y is 0. Then if 

(1) F (0, y) = 0 for y in an open interval containing 0, and 0 £.(R n 
is an isolated stationary point of F, 

n | 

(2) F is analytic in x andy in a neighborhood of (0,0) in |R x R , 

(3) A ( y ) = F (0, y ) has a pair of complex conjugate eigenvalues 
X and X such that 

X (y ) =a(y) + iw(y), (A.2) 

where a) (0) = u) Q >0, a (0) = 0, a' (0) # 0, (A.3) 

(4) The remaining n-2 eigenvalues of A (0) have strictly negative real 
parts, 

then the system (A.I) has a family of periodic solutions : there is an e H >0 
and an analytic function 
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y £ ) - £°° yH gi 

2 1 (0< e < e H ) (A.4) 

such that for each e £ (o. e ) th^r* ev ic+,- ^ • .. , . . ' 

’ H' ™ ere exists a periodic solution p £ (t) occurring 

H H 

for y - y ( e ). If y ( £ ) is not identically zero, the first non-vanishing coeffi- 
cient y . has an even subscript, and there is an ^£(0, e H ] such that y H ( £ ) 
is either strictly positive or strictly negative for e£ (0, £j). For each L > 2 it/ to 
J there is a neighborhood n of x = 0 and an open interval I containing 0 such 
that for any y £. 1 the only nonconstant periodic solutions of (A.l) with periods 
less than L which lie in n are members of the family p £ (t) for values of £ 

satisfying y ( e ) = y, e£(0, e^) . The period T^( e ) of p £ (t) is an analytic 
function 

T H (e) = (2 tt/w o ) [l + z“t H . e 1 ] (0<e<e H ) (A.5) 

Exactly two of the Floquet exponents of p £ (t) approach 0 as e 0. One is 
0 for c £(0, e^), and the other is an analytic function 

e H ( e) = 2” e” e' 1 (A.6) 

The periodic solution p (t) is orbitally asymptotically stable with asymptotic 

C 

phase if £}*"*( e ) <0 but is unstable if $^(e ) > 0. 

4.2 Numerical Techniques Used 

Numerical experiments on a dynamical system may be used to inves- 
tigate its global behavior and to locate limit cycles, if any. This can be done 
by numerically simulating the dynamical system on a digital computer. This 
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differential equations over a large number of steps. 

Standard library subroutines such as DVERK (in IMSL library) and 
D02BAF and D02EAF (in NAG library) are available for solving a system of 
ordinary differential equations. In general, Runge-Kutta methods, such as Verner's 
method (DVERK) and Person's method (D02BAF) can be used without problem. 
However, for stiff systems, these methods fail, and methods specially suited 
for stiff systems, such as Gear's method (a predictor-corrector method) have 
to be used. The subroutines D02EAF and D02EBF in NAG library' are based 
on this method and were used for integrating the systems based on models 


with delayed neutrons. 
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APPENDIX - B 

AH IMTRRACTIVe'pRDSHA^FOR COMPUTING 
'''HE LARGEST Lupumny exponent 


This la an interactive program 
Lyapunov exponent CELAMDA) of a 
of 41 ne ns ion U . 


for computing the largest 
continuous dynamical system 


The e in at ions defining a dynamical system and its variational 
equations are Integrated simultaneously over a given steo sl*e 
starting tfl.th given initial conditions and an arbitrarily chosen 
variation vector of unit magnitude. After a step of Integration, 
the variation vector is norma lizegto unity and the process is 
repeated. Tpe largest Lyapunov exponent is computed by calculating 
the ratio of the logarithm £to the oasft 2) of the ratio of the 
magnitude of the variation vector after and before integration, 
averaged over a large number of steps,. ^ ^ 

PROGRAM LYAP 
PARAMETER N*1,MN=6,NW=9 

Ms# of equations in the system t , _ 

Nn s 2*N s. total ♦ of equations to be integrated, 

EXTERNAL LOREN£,ROSLER,REACrR 

Subroutine containing the equations to b® integrated is to be 
declared in the above EXTERNAL statement. 

T integer meq,tw,ifail,i,j,k,io,mskip,nstep 

m step = n of stens of integration ,• 

M SKIP = # of stens to be shipped initially ,* 

ID s input /output rate of steps. 

REAL X, XRNn,TOL, STEP f ZN0RMl,2N0RM2» TERM, SUM, ELAMDA 

X = the independent variable ,* 

ZNORMi = norm of the position vector ; 

7, n h R M 2 * norm of the variation vector i 

ELAHhA = estimate of the largest LYaognov exponent , 

REAL f CNN) ,W(NM,NW) ,CCC24) 



r 

n 

f 


r 

r 

r. 

C 


oa^e 7 

v ~ arrav containmy the fi dependent var 1 ables and 
t heir v variations . 

r c * l n , u , c , d , e , f 
r D’' mo. vi a,b,c,d,e,f 


Values of the parameters 


,:.u,Q. 5 , 8 «' 0 e -S*V 


n ATA C»DrF,ij/*:,u,u,o#o»u,’ . 

WA C,<>,F,G/2,u,l ,6,i.0,fl,-2.O/ 

TYPE b 

POPIMflH , 'PI 
ACCEPT*. A 


nive A") 


r 

c, 

c 

5 


FsH-A 

npn^CUFXT=:23,DFVTCF='DSK 1 ' ,FlLE-'i-XAp,DAT'} 

Initial rondi. tions for nonlinear system 


TYPE 5 

for ha T c JH , 'Pi, aive toe initial conditions ') 

ACCEPT*, ( Y ( T ) , 1-1 , N) 

C.,, initial conditions for linear system (arbitrarily chosen 
C.„, variation vector of unit norm) 

C 

V C 4 )ss 0.75 
■ Y f i ail'' ' A . 

YC 6 )*s 5 RTC 3 ,/i 6 .) 

C... Integration tolerance, time per sten,t of integration steps, 
C.,,and t/n rate 
c 

1 JonHAIUH .'PI. C.IV* tOfc.STEP »STEP,la,NSKIP') 

ACCEPT*, TOL, STEP, NSTEP, ID, WSKIP 

r 

r . . .Tni tial ization for integrator 

r* 

MEOsNni? IWsMW „ _ 

X = 0.0,*XX==X+FLOAT(NSKIP) *STFP 

s U M 

TNhEX=0 . 

100 no 200 I=J jfbSTFP 

index=index+i 

TNDsi 

xend»x+step 

C, , .Integrate the DDE's by an IMSL routine usinq Verner's method 
^ CALL DVEPKCfiFQ,RFACrR,X,y,XEND,TOL,XHD,CC,NF^, W ,IEP) 



page 


2 

r 

c, 


cn Dl .IC^^ D * f JE- :i * OR »l F:R * NE *' 0 J I'YPe 2,Tf»D, IER 
FORhAlfiH .'Error checks IND S ,ii 'iers' , 13 ) 

.Calculate the norms of the position and variation vectors 


ZNHRKlsf ,o?ZKQPM?=o.n 


n D 2 C Jsl.N 


7 A 


ZNOK^isZt'ORMI +Y( J)**2 
ZNnRM2=ZNORM2+Y(N+J)**2 
CONTINUE 

7.HmR«isSORT(7NORMl) 

7iV0R^2 = S0RT£7hiORY2) 


Normalize the variation vector 


30 


DO 30 KsN+l.NN 

YCKJsyf K)/ZNhRM2 

CONTINUE 

TFClNDEX.hr. NgKIP) GOTO 50 


• • • 


calculate the largest (positive) Lyapunov exponent 


TERH*ALOC(ZNORM2)/ a {jDGC 2* ) 

SUMsSUM+TE'RM 

ELAHDAbSUM/CX*XX) 

c 

C,. .Print every 10 iterations 

SO IFCHOO(INOEX,IO).NE.O) SOTO 200 

WRTTE( 23 , 3 ) Xj.ELAMOA.ZNQP*! .ZNOR« 2 ,Y 

TYPE *,X,EhAM6A 
FORMATCIH , 10( 2X, Ell . 4) ) 

r* 1 T ■ , ■ . . . . . 

TYPE 3 .X,ELAMDA,ZN 0 RM 1 .ZNQRmZ, Y 

T i P E 4 

FORHAKiH 'How many more steps *') 

ACCEPT*. NSTEP 
IFC'JSTEP.GT.O) GOTO 100 


3 

200 

C 


END 


r+*****t* ****************************************************** *******? 

SUBROUTINE LORENZCN. X. Y.YDOD ! 

C ’ I 

C The Lorenz system (chaotic) to test tne Program . > 

REAL X,YC6).YD0TC6) | 

?8Rm wmiu »**rco.tc*» 

YUOr(3)BYCl)*YC2)-4.*Y(3) 

T r83 , r{8!:{teil;«*H4)-T(*i.m».t(„ 



':4 0 


^”” t *s”^ 5 : f^j ,, ^^|jj^**j***;j*’***»*‘********»**«»**»******** 

C The dossier system (chaotic) to test t*e program . 


PARAMETER A= 0 .t 5 ,B= 0 , 2 »C=in.h 

real K.x (6),YdnrU) 
Yoora 5 =“Yc^)-Yc 3 ) 
YDnr( 2 )=Y(l)+A*YC 2 ) 
Y00r(3)=B+Yp)*CY(l)-C) 

Y 00 r( 4 )=-YCS)-Y( 6 ) 

Y 00 f{ 5 )=Y( 4 )+A*Y( 5 ) 

Y00r(&)=Y(3)*YC4)+CYfl)-C)*YC6) 

R & T U R N 

e-mo 

SUBROUTINE RFACTR(N,X, Y f YDOD 


C 

r 

c 


c 

r 


Two temperature feedback model without delayed nutrons 
integer n 

REAL X,YC3).YDOTC3) 

REAL A,B,C,O r E,F 

common a;b,c,d;e,f 

The dynamical system 

YpOr^|)aAfYi[ 2 )+EfY( 35 tA*YCt )^YC 2 )TE*Y£l)*yf 3 ) 


Y 0 Ore 5 jaCC-O)*Y(i) 
YUOr( 3 )=F»YC 2 )-F*y 


„*YC2HD*Y(3) 

C3> 


The variational equations 

Y0Or(41 = CA*YC2)+E^Y(3))l'Y(4)fCA + A*YCJ))*Y(5) + (E+B*YCi)3YY(6) 
YD0rC5) = (C-DmC4)-C*YC5U0*Y(6) 

YUOrC6)=F*Y(5)-F*YC6) 

RETURN 

END 


Eigenvalues at a = 0 : 0, - c 



Fig. 1. Movement of eigenvalues in one- temperature 
model without delayed neutrons as a 
decreases (c>0). 



Fig. 3. Hopf bifurcation created in one -temperature 
mode.l without delayed neutrons by an extra 
parameter . 


2.0 



Fig. 2. Orbits for One- Temperature Feedback Model 
without Delayed Neutrons. 


Parameter Values 



Fig. 4. Orbits for One-Temperature Feedback Model 
without Delayed Neutrons and with Constant 
Power Removal. 



Notations : 

o Eigenvalues at a = 0 
***** Prohibited regions 



(a) Case a<0, c>1 



(b) Case a<0, 0<c<1 

Fig. 5. Movement of eigenvalues in one- temperature model 
with delayed neutrons as a decreases. 
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Fig. 6. Projection of Orbits for One-Temperature 
Feedback Model with Delayed Neutrons. 




(a) Case a<-a* (b) Case a = -a* 



(e) Case a>f 


Fig. 7. Movement of eigenvalues in two-temperature model 
without delayed neutrons for the case a+q = 0. 



Fig. 8. HOPF bifurcation in two-temperature model 
without delayed neutrons as the parameter 
a increases past the critical value a*. 



Parameter Values 


a = 4.5 j cj =—9.5 ,k = 2.0, 
p = 0.25, r = 1.0 



Fig. 9. Projection of an Orbit for Two-Temperature 
Feedback Model without Delayed Neutrons. 



Fig. 10. Projection of Another Orbit for Two-Temperature 
Feedback Model without Delayed Neutrons. 


TABLE 1 


TYPICAL VALUES OF PARAMETERS 


PARAMETER 

VALUE 


THERMAL REACTOR 

F AST REACTOR 

3 

0.0066 

0.0035 

X 

0.077 s' 1 

0.065 s _I 

1 

4.0 x 10' 5 s 

5.0 x 10~ 7 s 

P 

o 

1040 MW 

3000 MW 

a f ’ a c 

-10" 5 to + 10 -5 

-10" 5 to + 10~ 5 

K 

o 

1 

X 

c_» 

1 

10" 7 KJ' 1 

7 

1343 K 

1783 K 

T° 

c 

500 K 

600 K 

T. 

600 K 

700 K 

l 




In models without delayed neutrons, we take 1 to be of the order of the effective lifetime, 
i.e., 


1 


37 X 




“7-ft/U£CTg>/?/ES CALCt Pofr 
QbJE -T£MP££dTUftE MobBL WITH bEJ-AVEb M£UT*0*tS 

TOL= 0. 1 0E-10 2 
INITIAL COMO-ITID^S: 

O.SOOOEi+OO 0.500CE+03 0.2000E+00 


X 

XU) 

: .OCOOE+CC 

9.60395+90 

*.lOOOE-01 

D.29345+D0 

„ .20095-01 

9.90135+09 

V.3000F.-01 

0.4 9 6 7 5+ 93 

; .40005-01 

?. 49545+03 

..50005-01 

0.494 15 + 3 9 

v .60005-01 

9.4919503 

t .70005-01 

9.43935+00 

V .8000E-'01 

0 .'18 7 55 + 99 

, .93D0F--01 

0.43565+90 

’ .1 O0OE+OO 

0.4 3 355+99 

’ .lOOOEtOO 

0.43355+03 

. .-2000E+00 

0.46525+03 

L.30OOE+O0 

0 . 4 4 9 9 E + 0 9 

*.40005+00 

3. 4 -36 35 + 00 

r .5O0OE+OO 

9.42575*30 

1.60905+00 

0.41595+03 

l .7O0OE+OO 

0.4-97 35*03 

I .eoooE+oo 

0.99535-01 

: .9300E+00 

0.92 51 S'-SD 1 

■ .1 O0OE+01 

9.35335-31 

I .1 0005+01 

0.85335-31 

1 .2C-00E+01 

0 .4 3 55 5-0 1 

; .3C905+G1 

0.-22 4 35-31 

. .40905+01 

0.4 1535-01 

1 .59005+01 

0.60755-32 

; .60005+91 

0.31545-32 

; .70005+91 

0.45335-32 

r . boooe+oi 

D. 92 425-3 3 

: .90005+01 

0.51925-33 

.4000 5+ 0 2 

9. 27 2 05-03 

l .10005+02 

0.-2 7205 -3 3 

. .20005+32 

0.-71 475-35 

‘ .39005+92 

0.4 3 5 IE-35 

=.40995+02 

0.631 95-36 

..50005+02 

0 .64 3 1 5 -3 5 

r .60035+9? 

0.65255-3? 

r .79005+02 

0.4 0345-37 

l .PoO'jE + 02 

0.4 5 4 35-37 

, .9G9O5+02 

0.-25975-38 

*.10905+03 

0.32575-39 


Y(2) Y( 3) 


0.-50005 + 00 0.2000 5'+ 00 

D.-4972E + 00 0.-20025 + 09 

9.4912E+0U 0 . 2 0 0 3 5+ 0 9 

9.4913E+PC 0.2002E+03 

0.46R4F+PO 0.20025+39 

9.4855E+0P 0.2001E+09 

3.48265+00 0.2000E+09 

0.4797E+00 0.49Q8E+09 

3.4756E+0C 0 .4 997 E'+ 3 9 

0.4739E+OO 0.4994E+09 

3.4710E + 00 0.49925+09 

0. 47105+00 0.4992E+00 

0.4428E+OO 0.4953E+09 

9.4156E+0G 0 . 4 8 9 5 E+ 0 9 

9.3897E+P0 0.4624E+09 

0.3650E+OG 0.4 7 4 4 E + 0 9 

0.3417E + P0 0.4 659E+03 

0.91985+00 0.4 5735+09 

3.29915+00 0.-1 48BE+O0 

9.2798E+C0 0.4 4O5E+09 

0.-2616E+GG D.4 3265+09 

9.261 GF+O0 0. 4 3265+09 

3.1334F+00 0.-7 05 9 5-01 

3. -SR 315-31 0.36615-01 

5.3536E-C1 0. 4922 E- 01 

3. 16495-01 0.4 0055-01 

3. 95735-02 0.624151-02 

9.5107E-02 0. 2800 E- 02 

3.-2798E-02 0.45345-02 

3.1 5455-02 0.64775-03 

0.6236E-53 0.45225-03 

9 . 32 36E-C- 3 0.45225-03 

9.21615-04 0.4 1 8 B 5- 0 4 

C.40P3E-05 0.22465-05 

0.2514E-P5 0.43R35-05 

3. 94425-06 0. 61995-06 

3.1 959F-06 0.4 07 6 E-0S 

9.1235E-C6 0.6b62E-07 

3.46G1E-07 C. 25465-07 

3.7S98E-0B 0.42745-08 

9.'995?E-09 0.6487 5-0 9 



*TABL£2. (CCA/TD.) 

rOL= 0.103-02 
INITIAL COSKTIDn’S'I- 

0.15003>0l 3.<100DE+31 0.500OE+00 


X 

1(1) 

1(2) 

Y( 3 ) 

I .O0OOE+OC 

0.1530 3+ 0 1 

0.1 CD0E+01 

P. 60003+00 

;■ .-loooE-oi 

0 .*24 1 43 + DG 

3.9926E+00 

0.49633+03 

•: .2030E-0! 

0.24 3 7 3+ 03 

0.9853E+OO 

0.4 971 3+0 3 

1 .30003-01 

3. 2 '4 3 : 3'+ 3 : 

3 . 9 7 7 9 E + 0 C 

0 .4 87 9 3'+ 0 3 

r .4oOOE-tH 

0.23933+03 

0 . 9 7 0 5 F + D 0 

0. <4 837 3+ 00 

. .50305-01 

5 .'23 ? 5 3 + 3 3 

3.9&33E+0'' 

n. 47073+00 

. .600CE-01 

0. *23 793 + 03 

5.9551F+(H 

0. <47563+00 

: .70003-01 

0.23723+33 

:.3.?r9F40o 

n. <471 7 3 + 0 0 

. .eoooE-oi 

3.23543+03 

0.941 B F + 0 0 

0,4678 3 + O'* 

t .90033-01 
.lOOOE+OC 

0 .'?3 5 6 3 + 0 5 

0.'3343F+O0 

0.46393+00 

3.23193+0? 

5.9279E+00 

0.4501 3+00 


L.1000F+00 

0. <23195 + 00 

0.9279E+00 

0.-4 601 3+0 5 

t .7o30F+00 

O. 22 543 + 05 

5.661 5E+0P 

0. -42463+00 

.30003+00 

0. <71593 + 00 

0. 93953+00 

0.99323+00 

-.43003400 

0.20733+00 

3.74453+00 

0. <36523+00 

t .50003400 

0. <19573 + 00 

0.59793+00 

0.34053+00 

: .60003+00 

0. *13533 + 00 

3. 64533+00 

0. *31733+00 

: .700QT+00 

3. *1 7503 + 0: 

0.551 3F+ 00 

0.79653+00 

: .eoocF+cc 

3. <1651 E + CO 

5. *56353 + 03 

0.77733+O0 

*. .900GF+0C 

0.15553+00 

0.52253+00 

0.75953+30 ; 

.100CE401 

0.147 15+00 

?. 45743 + 00 

0.7 4 3 0 3’+ 0 0 

V .lOOOE+Ol 

3.44743+00 

0.48745+00 

0.74303+00 

E .2000E+01 

3 . <7 3 033 -I? 1 

0. *24533+00 

0.4 2723-+03 1 

.‘.30OOE+O1 

0.4355E-D1 

0.125? 5+ 0 5 

0.66523-01 

r.4000F401 

0 . <2 1 0 5 3-30 1 

0.64223-01 

0.-3466 3- TV 1 i 

.’.5000E401 

0.10333-01 

0.-3310E-01 

0.4 «0l3-01 | 

.600CE+01 

0.56733-02 

0.1721E-01 

0.93993-02 I 

■„ .70003 + 01 

5.3 0 07 5 -D 2 

0.91033-02 

0.49823-02 ! 

L .eoooE+oi 

0.4 5 S 55 -S 2 

0.-481 OE-02 

0.75373-02 I 

: .9000E+01 

0. <3*333-33 

0. *2 559 E- 02 

0.4 1043-02 i 

”.10003+02 

0.4 657 3 -.03 

0.14! 2E-02 

0.7 7 48 3!— 03 j 

: .10003+02 

0. <15573 -3 3 

0.141 23-0? 

0.77-183-33 

1 .20OOE+D2 

0.9 7 57E-I ; 5 

0.29553-04 

0 .4 62 4 3— 0 4 j 

I .30003+02 

3.-271 95-35 

0.62253-05 

0.45223-05 | 

0 D ^ p* ■f' 02 

0. <37 4 43-1! 5 

0.11 253-05 

0.62063-05 | 

'I l5O0OE+h2 

0.4 5973-^5 

0.4? 983-06 

0.76453—05 j 

. .60003+02 

0.< 25 4 2 5-1*7 

0. 75993-07 

0.42023-0/ ! 

.70003+02 

0.91 OlE-23 

0.2 4 2 23-07 

0.43403-07 

r .eoooE+o? 

0.19713-03 

0. 5901E-08 

0. <32663-03 

' .9000E+02 

0. 10 5 3503 

0. *31553-08 

0 . <1 7 4 6 3- 3 3 

.10005+03 

0. 13453-09 

0. -40913-09 

0.72523-39 



'TABLE. 2 . (COHJJ) 0 


TOL= D.10E-102 

initial pbmross: 

-0. '8 000 5'+ 0 0 -”>.<5D00E+30 D.5000E+00 


X 


rci) 


X(2) 


y(3) 


; . 0 a 0 0 E t 3 D 
-Vl C0OE-91 
! .zopoe-oi 
• .3ODDE-01 
t'.4oOOE-:01 
:‘.50!>OE-Ol 
. .6O00E-01 
.7O0OE-01 
. .800OE-O1 
. .9O00E-O1 
.lDOCE+OP 


; .1 OO0E+O0 
\2000EfD0 
k .3O0OE+O0 
.4000E+ 00 
r .5000E+0O 
.60005+ DO 
...7O0OE + 00 
: . 8 0 DOE + oo 
;,.90 00E+00 
.1 OOOE+01 


- .-1 3O0E+O1 
.2&OOE+01 
. .3O0OE+O1 
„;.4O0OE + D1 

* .5Q0CE+O1 
: .60005+01 
i .70OOE+O1 
: .BOOOEi-Ol 

* .00OOE+O1 
‘*.1 O0OE+O2 


: . 1 00OE+ 02 
. 2 0 D 0 E + Q 2 
: .30005+02 
; .40005+02 
! .50005+02 
.60005+02 
' .70005+02 
" .60005+02 
; .90005+02 
' .10005+03 


'3 • <3 3 O 3 5 + 0 3 
0. 63505+30 
■0 . 5 3 2 4 E‘+ 0 [.» 
■•0.6733E + OO 
■'0.67 5 15 + DO 
0 .67 1 45 + 00 
•■0.6575 5 + 0 0 
■‘0 .6537 E + OD 
■•0.65995 + 00 
•’0.65595 + 03 
•0.65 2 05 + 00 

■0.65205+00 
•’0.6375 £ + 0 0 
•0.65255+00 
■0 . <5 1 2 0 £ + 0 0 
•3. <15955 + 00 
••0.43 7 05 + 0 0 
•3.65 7 3 5 + 0 0 
■O.'Bl 1 35 + 0 3 

• 0.<271 95 + 30 

• 0. <23795 + 00 


-0.50005+00 
-P.501 9E+0C 
-,0 . 5 0 3 7 F + O 0 
-D.5055E+00 
-P.5072E+00 
-P.5D38E+0C 
-D. 51 045+00 
-D. 51205+00 
-P.51 35E+00 
-0.51 49E+0C 
-0.5153E+OO 

-D. 51535+00 
-P.5272E+00 
-D.5328E+O0 
-0.5332E+OO 
-0.5237E+OC 
-P .51 96E+00 
-P.50S5E+00 
-0.4901 E + OC; 
-0.471 25+00 
-0.4506E+0C 


3 *<23 7 95 + 0 0 
■3. <33735-01 
• 0.43 3 45-01 
•3.-221 95-91 
•9. <11 3? £-01 
•9.5333E-02 
•9. <3 3 9 7.5 -02 
•3. <15275-92 
•3 .<3 7 2 35-0 3 
•3 .47 045-33 


.3. <47 045-03 
.3.97525-35 
.9. <271 95-35 
•3 .37 055-9 5 
.3.45295-05 
.3.-2 4375-37 
.9.71155-33 
-3.<l 91 05-0 8 
-3. *19 2 2 5-3 8 
.3.1-35 05-39 


D.-4596E+00 
0, 24935+00 
•0.1 2B9E+QO 
•0.6654E-O1 
•D.<3440E-01 
■0.1786E-O1 
■0. *93345-02 
•0.49255-02 
•P.2645E-02 
•0.14275-02 


•0.1427E-O2 
•0. <29535-04 
0.82555-05 
■0.1 131E-05 
-P.4572E-06 
•D.7439E-07 
■0 .<2l 96E-07 
■D .5 8 5 95 -08 
•0. <31 335-08 
■D.4075E-09 


0.50005+00 
0 . 4 B O 2 5 + 0 0 
0.45095+00 
0.441 95+00 
O.42335+00 
O.4O5O5+O0 
0.-3 871 5+0 0 
0. <3 69 6 5+ 0 0 
0. <35245+00 
0. <33565+00 
0. <31915+00 


0.< 3191 5 + 00 
0.<1 7245 + 00 
0.65315-01 
-0. <3 60 2 5-01 
-0.1 0535+00 
-0.1 557 5+ 00 
-0.4 9055+00 
-0.21255+00 
-0.22455+03 
-''.<22835 + 00 


-0.22895+00 
-0.1 4455+00 
-0.73785-01 
-0.-37 1 9 5-D t 
-0.1 9O95'-0 1 
-0.0874 EV 0 2 
-0.61745-02 
-0.27115-02 
-0.1 4565-02 
-0. <784751-03 


-0. <78475-03 
-O.1627 5-04 
-0.4 531 5!-0 5 
-0.61955-05 
-0.25575-05 
-0.40865-07 
-0.1 1925-07 
- 0 . <3 1 9 4 5!- 0 3 
-0.1 709;— 03 
-0.2235 5- O 9 



-TABl-E Z.CCoNTb.) 


rOi,= 0.10E-D2 
initial comoi riass: 

-0.60005*03 -:3.633OE+30 -0.6OCOE+OO 


r(i) 


r(2) 


1(3) 


' .O'OOOE* 0O 
.■I 0005-02 
■* » 2 C* 3 O £ — 0 2 
. .30 00E-02 
, ,4'iOOE-02 

• .5? 3 0 2-02 
.60002-02 

. .70305-02 
.gO 305-32 
.90305-02 

.10332-01 
.23305-91 
.33 30 £- 31 
. 4 3 0 DE-.J1 

* .5033E-01 
; ,6900E-:31 
i .7030E-:01 

« .eoooE-oi 
v .93 DOE-31 
...130OE+OO 

3.10335+00 
I .23305+00 
.3033E+00 
.43005+03 
.53035+00 
.63305+30 
.73035+00 
.€3005+30 
.93035+00 
, .10305+01 

4.10305+01 
..20 DO 5 + 31 
.30335+31 
: .4030E+31 
-.53305+31 
.63335+01 
.79005+01 
, .8030E+01 
, .90305+01 
.1 303E+02 

.10335+32 
.■23305 + 32 
.30305+02 
.40305+32 
" .50305+32 
, .63305+32 
, .70335+02 
. .90DOE+32 
.90335+32 
.13305+03 


O .'5 333 
•P.'l 233 
3.1433 
3. '2933 
3 . '3 7 3 3 
3.4 1 53 
3.4371 
3.4155 
3.4455 
3.4134 


£ + 33 
£ + 3 3 
5 + 3 3 
5+33 
5 + 03 
5 + 33 
5 + 3 3 
5 + 33 
5 + 3 3 
£ + 3 3 


3.43325+33 
3. '377 95 + 3 3 
3. '32 4 35 + 0 3 
j. *23135 + 33 
3 . 21575+03 
3. 2 '1182 + 33 
3 .1 33 3 5 + 0 0- 
3.1 51 75+ 33 
■3.14 415 + 0 0 
3.12595+03 

3.12595+03 
3.1652 5-:3 1 
-:3. '31 9 55.-31 
-;3.o5135-:31 
-;3. <591 3 5-31 
-p.75 3 12-01 
-3.7 7 155-31 
-•3.7 7155-31 
-;3 .7 5 3 25-31 
-;3. '733 42-3)1 

-:3.'7 3 31 5-01 
-9.41395-01 
-9. 21 2 25-31 
-9.139 32-91 
-;3.'55 4 3 5-32 
-9.29 5 35-9 2 
-9.153 25-3 2 
-9. '3 13 35.-33 
-9.472 42-93 
-3. 25 3 4 2.-9 3 

-0.25 3 4 5-03 
-9.73912-3 5 
-9.11355-25 
-9.7 7 132-2 5 
- 9 . *36 2 12.-3 5 
-9 .'55 3 32.-3 7 
-9. '37 7 5 5-37 
-9.29 3 ! 2 -9 7 
-9.253 3 2-2 3 
-9 .'34 2 4 5-3 9 


-D.-5033E+00 
-3.59975+00 
-9.5991E+03 
-D.5933E+09 
-0. 59715+00 
-D.5934E+00 
-D. 59535+05 
-9.59432+03 
-9.5933E+30 
-P.5922E+30 

-D .59 1 25+00 
-P.581 4E+ 00 
-0.5721 E+OO 
-D.5634E + 00 
-D.5552E+0O 
-0.5473E+OO 
-P.5399E+00 
-D.5327E+00 
-D.5259E+0C 
-0.51 93E+0C 

-0.51 93E+00 
-D.4610E+00 
-D.4239E+00 
-D.'3853E + 00 
-D.'3518E + 00 
-D.9230E+00 
-0*’3O4iE+OG 
-0.2326E+OG 
-0.263OE+OD 
-0.2451 E+ 00 

-0.24 5 IE* OO 
-0 .1 2 16E+O0 
-P.64D2E-01 
-0 .'3296E-01 
-0.1739E-O1 
-0.-3991E-O2 
-0.4832E-C2 
-0.2552E-O2 
-D.1427E-02 
-0.7633E-O3 

-P.7683E-03 
-D.2243E-04 
-P .'36 3 lE-05 
-D.2353E-05 
-0.11 02E-05 
-D.'17D8E-06 
-0.11 52E-06 
-D.-5151E-07 
-p .7 9 l 3E-3 8 
“D.1031E-08 


0.6000 

0.5996 

0.5985 

0.5971 

0.6955 

0.5939 

0.5922 

0.5905 

0.5887 

0.5370 

3.5853 

0.5690 

0.5537 

0.5394 

0.6261 

0.6135 

0.6015 

0.4902 

0.4795 

0.4592 

0.4692 
0 . *3 8 7 2 
0.-3292 
0.2356 
. 2 5 1 5 
0.2240 
0 . 201 3 
0.1 323 
0.1 661 
0.1 520 

0.1520 
0 . *7 1 4 3 
0.35B6 
0.1829 
0.9443 
0.4958 
0.2644 
0.1 404 
0.7848 
0.4222 


540 3 
5*0 0 
5 + 0 0 
5*0 0 
54-0 0 
5*0 0 
5 + 0 0 
5*0 0 
54 0 0 
5*0 0 

5*0 3 
5*0 0 

£:+oo 

540 0 
54-0 0 
5*0 0 
54-0 0 
£400 
5400 
54 00 

140 0 
:4oo 
:400 
:4 03 
:40 0 
•4 0 0 
:4oo 
:4 0 0 
:4O0 

:*oo 

5400 
5-01 

5401 
5401 
5-02 
5432 
£402 
£402 
5403 
5403 


■0.42225403 
•0.12325404 
■0.1 97 7 5-0 5 
■0.1 2915405 
•0.6040 5405 
*0.93205407 
•0.6341 5407 
•3. '3 36 2 540 7 
■O.42875-03 
■0.564 35409 



"TABLE 2 . (cohJ'Tb.') 


rDL= 0.10E-D2 

iNiriftb cooiirois: 

9.60002333 -P.oOOOE+OO -0.4000E+00 


XC1) 


Jf{ 2) 


XC3) 


.00002+00 

.1D30E-02 

. ,2:dde-:o2 

.300PE-0? 

. .4000E-02 
. .503 3E-102 
.60 302 — 02 
.7 0 O0E-.O2 
. 8 C 3 G E-.3 2 
_ .9030 E- 02 

, .1C OOE-Ol 

: .scooe-ioi 

. . 3030E-01 
i- ,4oooE-:ai 
I .5O0OE-101 
P.6OO0E-1O1 
K700DE-01 
s .8000E-01 
t .9000E-101 
.10OOE+OO 

. .40302*00 
.2l'DOE*00 
v. .SdDOE + CG 
.400OE+C? 

, .SDCOEtOr- 
. .eeooE+oa 
. 7 0 0 0 E ♦ 0 0 
i.8QOOE+OQ 
.’-.•9000E + 00 
: .lOOOE+Ol 

.ICOOE+Ol 

. .20 COE* 01 
■ •.3C30E+01 
, .4fODE*01 
.5 0 3 0E + 01 
.6COOE+01 
, .73002+01 
.6D&DE+01 
l-.9C0O2+'01 
-.4 COCE+02 

. » 1 0 0 0 E + 0 2 
.20002+02 
.. .3GOOE+&2 
.40 0CE+ 02 
.50002+02 
. .60 OC-2 + 02 
■ .70505*02 
. .6D33E+02 
-.9DO0E + 02 
; .1000E+03 


-P 

-P 


3.63032+30 
G.'i332E>03 
3. *13 41 2-01 
-P.-3-3 3 3 2-132 
-P.'!9l5 2-lJl 
■0 .-2 '4 5 5 2.-131 
-13. -25372-101 
-p. '27 9 3 2-01 
•p. ‘2 ‘3 4 55-01 
•p. <29 9 0 2-01 

-0 .'29 3 9 201 
•P.3414 2-B1 
O .'33 332-01 
•P. -4 2 7 32 -01 
'3.4555 2 >B 1 
*.•5014201 
.•5 3 45E01 
p. 56 5 32-101 
O.o 9 4 1201 
O. *>235201 

O.52 05EO1 
0.3 3 51 EOl 
O .3913 El-D 1 
0.32 22E01 
-0 ,9 2 2 7 El-30 1 . 
•P. '9 3 362.-01 
0.37 41201 
0.3371201 
O. '7959201 
0.-7 5 5 3201 

O .'75532-101 
0.403 7201 
O. '2055201 
0.4-952201 
O. *5 4 93202 
-0. <28 3 120 2 
0.4535202 
0.3 3 5 7203 
0.4 5 02 20 3 
-p. <2 '3 6 1 2-3 3 

0.2 '3 5 1203 
-10.66 1 2 205 
O .4 9 7 3 SO 5 
-p .63 212.0*5 

O .14 3 5 206 
-13 .6926207 
0.9 7 312.0 3 
-S0.47 3 52-D3 
0.6 3 35209 
0.3 42 22-1 3 


-P.5030E400 
-P.-4991E + 0Q 
-P.4935E+00 
-P.4930E40C* 
-P.4975E+00 
-P .4971 E+ 00 
O.4956E+00 
-D.49S1E + 00 
-P.4957E+00 
0.4952E400 

-P.4947E+00 

-P.4901E+OC 

-P.4855E+0C 

O.4812E+0G 

-P.4758E+00 

-D.4725E+00 

-P.4S94E+0Q 

-D.4642E+00 

-P.4602E+00 

-P.4552E+D0 

-P . 4 5 5 2 E 4 0 0 
0.41 97E+OC 
-P . 3 8 7 8 E + 0 0 
-P .'35 97E+O0 
-P.-3343E + 00 
-P.31 12E+00 
-P.2901E+00 
-D.-2707E + 00 
-P.*2528E + 00 
-P.-2351E+O0 

-P.-2361E+00 
-0.121OE+OO 
-P.5221E-01 
O.-3214E-01 
-P.1653E-01 
O.3721E-02 
-P.4651E-02 
O.2535E-02 
-P.4352E-02 
0.71 74E-03 

O.-7174E-03 
-P.-2007E-04 
-P.6335E-05 
O.2374E-05 
-P.4298E-06 
0.1816E06 
O.3022E-07 
-0.14 51E-07 
-P.1951E-08 
-P .'2539E-09 


•0. 400023-00 
*0. 398721+00 
•0.397923-00 
*0.3 972 21+00 
•0.3 965 2+0 3 
•0.39592+03 
•0.395321+00 
•0.394721+03 
•0.39412+00 
•0.39352+00 


•0.3929 

•0.3869 

•0.3810 

•0.3753 

■0.3698 

*0.3645 

*10.3592 

*0.3542 

-0.3492 

*0.3444 

•0.3444 
*0.3023 
*0. <2688 
-0.-2415 
•0.-2185 
*0.4 990 
-0.4 821 
•0.1673 
•0.4542 
•0.142 4 


21+0 0 
2+0 0 
21+0 3 
21+00 
21+0 0 
21+00 
21+00 
2i+00 
21+00 
21+00 

21+00 
21+0 0 
21 + 0 0 
21+0 3 
21+0 3 
21+0 3 
2+0 3 
2+0 3 
21+0 0 
2+00 


• 0.4 4242+00 

• 0.691121-01 

• 0.3 4812-01 

• 0.4 7832^-01 

• 0.-918 42-0 2 
•0.48072-02 
•0.-2565 2-02 
*0.1 3952^02 
-0. '748921-03 
*0.39432-03 

•0.394321-03 
■0.1 1022-04 
■0.32952^05 
■•0.11382^05 
-0. -23512^05 

• 0.-992521-07 

• 0.4 6452.-07 

• 0.<79552r03 
-0.4 06121-03 

• 0.4 3922-09 



-rABLB 3 . A -rRAJEc-TORy qA LCULA r&b FOR, 

-r\fjo- te-MPbRatur. b model without x>blavei> MEQtROM S 

4.500000 , 0.O500003E+00, 2.000000 , 0. 5000000 , 

TOL = 0.10E-D2 
INITIAL CDOirpNS: 

-0.«000i>00 *0 .'S OOOE+OO -0.8000E+00 


rtt) 


rt2) 


r(3) 


i. 1 0005+00 

".200OEt00 
„.3000~t00 
^.-lODOE+OCt 
. .5000E+00 
'■.6000E+0O 
, .70002+00 
.>€0005+00 
0.90002+90 
.10005+01 

, .1 1 005+01 
L. 12005 + 01 
- .1 3005+01 
14005+01 
>15005 + 01 
.*15005+01 
, .17005+01 
18005+01 
.19005+01 

i.2200E+01 
{. ;. 23005* 01 
...2430E + 01 
. . 250 0E* 01 

”1583?:?,' 

28005+01 
x ,. 29005 + 01 
.30005+01 

*.31 0 05 + 01 
, .32005+01 
.33005+01 
.34005+01 
.35005+01 
.35005+01 
.37005+01 
,.38005+01 
*. .39005 + 01 
-.40005+01 
: .41005 + 01 
42005 + 01 
.43005+01 
'.44035 + 01 
.45005+01 
^.4b005+ 01 
.47005+01 
.48005+01 
*.49005+01 
. .50005+01 
.51005+01 
: .52005 + 01 
^ 53005+01 

.54005+01 
' .55005 + 01 
..56005+01 
• .57005+01 
V .58005 + 01 
: .59005+01 
.60005+01 

■' .61005 + 01 
;. 62005+01 
.63005+01 
t .64005 + 01 
.65005+01 
.66005+01 
r. 67005 + 01 
*.66005+01 


-’0 . *7 0 1 5 E + 0 0 
-:3. <55355+03 
-O. 3321 5+00 
-0.55415-02 
0.45335+00 
' 0.4 1 045+01 
0.49095+01 
0.27735+01 
0.34355+01 
v-0. 35255 + 01 

0.23435+01 
0.46535+31 
0.64535+00 
-0.49355+00 
-■0.38935 + 00 
—0 • ( 7 8 2 5 5 + 0 0 
-O. <37505 + 30 
-•0.92 055 + 00 
-*0.94425 + 03 
-0.9557 5+00 
-0.95325+00 
-0.9651 5+00 
-0.95545+00 
-0.95455+00 
-0.95055+00 
-0.9537^+0: 
-O.9435S+00 
-0.92945+00 
-0.90535+00 
-0.97425+00 

-0.92725+00 
-O. ‘75325 + 03 
-O. <55705 + 00 
-0.90965+00 
-O. ‘29235 + 00 
0. 45715-01 
0. *44435 + 00 
0.40105+01 
0.46345+01 
0 .'2'4 02 5 + 01 
0 . ‘29 2 7 5 + 01 
u-D. <29992 + 01 
0. *2*4 7 5 5 + 01 
0.45445+01 
0.60005+00 
-•0.95215-D1 
-!0.'|9 3 05 + 30 
-]0. <71045 + 00 

:3:8iHU33 

-p. <9 14 05+00 
-0.93135+00 
-0.94045+00 
-0.94415+00 
-0.94395+00 
-O .94025+00 
-0.93235+00 
-0.92095+00 
-0.93315+00 
-0.97712+00 

-0.93955+00 
-O. ‘73525 + 00 
-O. <70535+00 
-0.69375+00 
-0.4312 5 + 0 0 
-O. '20025 + 00 

0.42145+00 
0.65375+00 ” 


-D.7935E+0O 
-D. 77135+00 
-P.<72755t00 
-P.65315+00 
-P.53555+00 
-D.35955+9C 
-P.10975+03 
0.21965+00 
0.60905+00 
0. 99705+00 

0.42392+01 
0. 14085+01 
3.13535+01 
0.11345+01 
0.95755+00 
0.74535+00 
0.5365E+O9 
0.3 47 85+00 
0.18005+56 
0.32285-01 
-P .<97 1 95-0 1 
-D.-21 035+C0 
-P.3D3BS+O0 
-D. 39435+00 
-P. 46335+00 
-D.6319E+00 
-P.58515+00 
-D. 63155+00 
-P. 66365+00 
-P. 69735+00 

-P. *71705+00 
-0.7255^+00 
-P.72355+00 
-0.-70415+00 
-P .66375+00 
-P .69355+00 
-P. 48365+00 
-P. <32185+00 
-P.97565-01 
0.19085+00 
0.52425+00 
0.3579E+00 
0.1 1042+01 
0.1218E+01 
0. 11875+01 
0.10525+01 
9.96745+00 
0.66945+00 
0.47 3 4E+ 00 
D.3028E+00 
0.14522+00 
0.6472E-D2 
-P.1175E+00 
-0 .‘22535 + OC 
-P.31875+00 
-P.3999E+00 
-P.459BE+00 
-P. 52925+00 
-P.67995+00 
-D. <61 935+00 

-P.65025+00 
-P.67105+00 
-P.6B035+00 
-D. 67 585+00 
-P.65365+00 
-D. 60325+00 
-P. <53225+00 
-D, 4i5?5+&0 


-0.7986 
-0.7397 
-0.7671 
-0.7239 
-0.6505 
-0*6354 
-0.3643 
-0.1219 
0.18*5 
0.541 3 

0.6898; 
0.1 t 5 5 = 
0.1284; 
0.1273 = 
0.1157; 
0.9833= 
0.7863; 
0.6890= 
O.<*0?9; 
0.2333= 
0.6146; 
-0*5287= 
-0.1709; 
-0.7740= 
-0.3638; 
-0.441 6= 
-0.5086; 
-0.5558; 
-0.6139; 
-0.6534= 


r+oo 
<*+0 0 
'+0O 
<*+00 
*+00 
r+oo 
:+oo 
;+od 
::+oo 
:*+oo 

:+oo 

:+oi 

;+oi 

;+oi 

:+oi 

:+oo 

:+oo 

:+oo 

:+oo 

:+0O 

:-oi 

r-01 

:+oo 

r+oo 

:+oo 

:+oo 

:+oo 

:+oo 

:+oo 

:+oo 


-0.68415+00 
-0.705/5+00 
-0.71695+00 
-0.71545+00 
-0.69785+00 
-0.65875+00 
-0.59065+00 
-0.43392+00 
-0.32775+00 
-0.1 1385+00 
0.15675+00 
0.46335+00 
0.761 35+00 
0. <9923 5 + 00 
0.11135+01 
0.11153+01 
0.10245+01 
O.67615+O0 
0.70225+00 
0.52345+00 
0.35205+00 
0.1 9385+00 

-0.1 3785+00 
-0.78535+03 
-0.37095+00 
-0.44435+00 
-0.60705+00 
-0.65965+00 


-O.6026 
-0.6358 
- 0.6586 
-0.6699 
-0.6b71 
-0.6468 
- D .633 4 
-0.6300 


5+00 
5+O0 
5+03 
5+00 
5+03 
5+93 
5+0 0 
5+00 ■ 


.500000 



-rft&LB 3. (c.e>fil-rC>.) 


;.* 900 Et 01 
, . 73005+01 
71005+01 
" . 72005+01 
{ . 73000+31 
,. 74005+01 
:,. 7500 E +01 
-^S&OOE + Ol 

2 . 77005+01 

k. ' 7800 E +01 
, . 79005+01 
; i; . 80005+01 

£ 1005 + 01 • 

*. 82005+01 
r‘.e 300 E+oi 
"'.« 400 E +01 
'•I. 85005+01 
„• . 86005+01 
. 87905+01 
" .£ 5005+01 
. 8900 E +01 
. 90305+01 

v . 91 005+01 
„ .« 200 E + 31 
. 93905+01 
, . 9400 E +91 
. . 95005+01 

l . 9 & 005+01 
'.. 97005+01 

. 95005+01 
,,. 99005+01 
. . 1 0005+02 
... 1010 E + 02 
1 . 10205+02 
. 1 0 3 0 E + 0 2 
Js.l 0 4 0 E + 02 
10505 + 02 
, . 1 060 E +02 
.1 070 E +02 
i-.'l OgOE +02 
f . 1 0905+02 

hn?si:8i 

~ .4 1 lCE+02 
r.l 1505+92 
j. .11 505+02 
.11705+02 
' .11805+02 
, .11905 + 02 
■ .12005+02 

i. 12105+02 

t . 1 220 E +02 

f . 12305+02 
. 12405+02 
f*. 12505+02 
*- . 1250 E +02 
12705+02 
280 E +02 
•, .1 250 E+ 02 
f .13005 + 02 
f; .131 0 E + 02 
; . 13205+02 
. . 13305+02 
u . 1 3 4 0 E + 0 2 
.. . 13505+02 
,. 13605+02 
. 13705+02 
». 13805+02 
T. 1390 E +02 
14005+02 

i, .141 0 E +02 
i’«l 4205 + 02 
. . 14305+02 
. .14405 + 02 

•' .14505 + 02 

j . 14605+02 
. 1470 E +02 

' . 14805+02 
. . 14905+02 
' . 15905+02 
. 151 05+02 
. . 15205+02 
. 15305+02 
’ .15405 + 02 
‘ .15505 + 02 


-0 .69 3 3 £*03 
-?3.*3 3 5 4£fD0 
-0*‘3537=:>OD 
-o.09S3£>ao 
•2.9137£>0D 
«P. *72252 + 00 
-O.7255E+03 
-0.72B5 £+00 
•’D.^17^s:>DS 
-•3.9053S + D3 
-’0.*3372£O3 
-:3.0S35 2 + 00 

-p.02?B2>00 
-P. 75722*00 
-0.6335203 
.3.67515+03 
-O.4l5S£+00 
-*0.4929200 
3.4 4445+33 
3.62235+30 
3.4 3312+01 
D. 46352+01 
9. *21 492 + 31 
w0. <2*4 3 7 2 + 31 
0.2‘4^ 3 2 + 31 
0.49545+31 
3. 4214 E> 01 
3. 4 ‘59 35 + 0 3 
•0 .9 3 552-2) 1 
•0 .45 1 15+03 
-D. 65332+00 
.0. *77455 + 33 
-0.03952+03 
-0.67525+33 
-0.69552+33 
-0. >93575 + 33 
-0.93992+03 
-0.90745+33 
-p.09942+00 
-0.08542+05 
-0.96332+03 
-0.93215+03 

-0 .7 8 572 + 33 
-0. ‘72135+00 
-0.62992+00 
-0.49995+00 
-0. *31 7 95 + 0 0 
-0.66945-01 
0. *25992 + 33 
0. *73215 + 00 
0.4^1 72 + 01 
0.47575+01 
0 .'22 312 + 31 
1-6. *23925 + 01 
0. *21732 + 01 
0.45935+01 
D. 97435+33 
3. *22755 + 03 
-0. *23 4 35 + CO 
-0.62235+00 
-0.69022+03 
-0 . */ 3 5 3 5> 0 3 
-O. 94372+OC 
-0.97175+03 
-0 .38352+03 
-0.99515+03 
-0.99555+33 
-0.99135+03 
-0.97912+33 
-0.35975+33 
-0.33395+03 
-0. *7 3 9 25 + 00 
-0. *72082 + 33 
-0 .6 4585 + 30 
-3.62772+03 
-0. *36235 + 30 
-’3.4 *3 592 + 3 0 


-P. *2 4952+30 
-E.-2M5E-01 
3.94942+00 
0.65362+00 
3.3 3 54E + 00 
3.43362+01 
0.4 1 ! 12+01 
0.10632+01 
0.93182+00 
0.-7592E+00 
0.67532+03 
0.39982+00 

3. <23732 + 00 

3.9043F-31 
-P.39862-01 
-r. 1 5452+30 
-10 .*25 5 02+ 0 C 
-D. 34222+00 
-P.41 75E+33 
-D .‘*81 72+00 
-D.33552+00 
-D. *§7952+00 

-13.61 342+00 
-D. 53572+00 
-P.64812+03 
-D. 64542+00 
-P.62532+00 
-D. 68202+00 
-P.5093E+03 
-D. 39862+00 
-P.24142+00 
-D. 32182-01 
3.22442+00 

3.50522+00 
3.76912+00 
3.9587E+00 
3 .1 036E + 0 3 
3.43312+01 
3. 38512+00 
3.7253E+Q0 
5 .65 342+00 
3.35442+03 
3. <22 7 12 + 00 
3. *84572-31 
-D. *42562-01 
- 0.4 5492+00 
-P. *25322 + 00 
-P. *33 8 62+0 0 
-P.4122E+00 
-D .47472+00 
-0 • 5 ? 5 6 E + 0 0 
-0.56832+00 

-0 . 5 9 9 4 2 + 0 C 
-0.61922+00 
-P.625OE+0O 
-0. 61732+00 
-P.5890E+00 
-0.53582+00 
-P.*4535E + 0 0 
-D. 32522 + 00 
-P. 45312 + 00 
3.67062-01 
3.32422+00 
0.58R7E+OC 
0.31432+00 
0.9533E+00 
3. *98382 + 00 
3.91792+00 
3.78S7E+00 
5.62972+00 
3 .46 43 E+OO 
3.30442+00 
3.45712+00 
0 **24392-01 
-0.94272-01 
-0.19852+03 
-0.28962+06 
-0.36832+00 
-0.4356E+OO 
-0.49232+00 
-D .52782 + 00 
-0.5^332+00 
-0 .59702+00 
- D . 53832+33 
-D. 53482 + 00 
-P.5&33E+00 
-0.5381 2+0-0 


-0.4178 
-Q.*2582 
-0.4667 
0.2108 
0.4902 
0.-7481 
0.9360 
3.4 022 
0.4 307 
0.-9127 
0.7711 
0.60R3 
0.4423 
0.-2B27 
0.4 355 
0 „ *2 4 Q 1 
-0.4 159 
-0 .*220 i 
-0.0114 
-3.0903 
-0.4579 
-0.61*39 


! ; :33 

2-31 
2+00 
2 + 30 
2+0 0 
2+0 0 
£+31 
£+01 
£’+ 0 0 
£+0 0 
E'+DO 
£ + 0 0 
£+0 0 
£+00 
£-0 2 
=*4 3 7 
£+0 0 
£+0 0 
£+0 0 
£+3 0 
£ '+ 0 0 


-0.5d1 52+03 

-0.697 9£+00 
-3.623 4 £+ 30 
-0.6369£'+ 0 0 
-0.6361 £+00 
-0.6l78£+00 
-0.6 770 £+00 
-0.6371 £+00 
-0.4007 £ '+ 00 
-0.i?501£+00 
-0.6194 £-01 
0 . 4 8 7 8 =1 + 0 0 
0 .4 472 £’+ 0 3 
0.6890 £+00 
0 .8 56 7 £ + 00 
3.9533 £+00 
0.9 46 5 £+00 
0.«543£+00 
6 . 7 3 4 5 £!+ 0 0 
0.681 7 £+00 
0 . 4 2 3 4 £'+ 0 0 
0.71P2E+O 
0.4 2 7 7 £“+ 0 0 
-3 .4 7 33 £-0 2 
-0.4 1 73 £+3 0 
-3 .9 1 °3 £‘+ C* 0 
-0.3 OP 5 £+3 0 
-0 .'3855 £+ 3 0 
-3.451 3 £+ 00 
-0.6062 £+00 


-0.6575 
-0.6842 
-0.6362 
-0.6152 
-0.6095 
-0.5925 
-0.531 8 
*0.4570 
- 0 . 0 2 ° 8 
- 0.1659 
0.4t34 
0.9802 
0.5233 
0.7315 
0.8683 
0.9149 
0.9783:: 
0*7001 = 
0.645b* 
0.4954= 
0.0437= 
0.4 988 = 
0.6505= 
-0.5579= 
-0.4633= 
-0 .*257 ? = 
-0.0402= 
-0.4108= 
* 0 . 4703 = 
- 0 . 51.89 = 
- 0 . 6567 ! 
- 0 . 5833 = 
- 0 . 6966 = 
- 0.69525 
- 0 . 5757 = 


£:+ 3 3 

£+33 
£+37 
£+33 
£+33 
£+0 3 

£*+55 

r+oo 

£'+ 0 j 
£-31 
£ + 03 
£+0 2 
=*'4 n o 

£’+33 
£!+ 0 3 
£’+0 3 
£ + 05 

£:+3 3 
£+03 
£+0 3 
£+03 
£-01 
£-01 
£+03 
2 + 0 3 
£+0 3 
£ + 0 3 
+ 33 
+ 0 j 
+ 03 
’+ 0 3 
: + 3 3 

: + o :• 
:o j 



-TAKL.S k.CtoyrrDd 


f«.77505t02 
f ..77 50E+ 02 
; .7770Et02 
■ >.77805+02 
, .37905+02 
,'•3® 705 + 02 
'■ .78105 + 02 
. .78205+02 
.*.78305 + 02 
! .73405+02 
*.76505 + 02 
J;. 76605+02 
. .76705+02 
".78805+02 
; .78905+02 
.79505+02 
i. 79105+02 
.79205+02 
.79305+02 
..79405+02 
.79505+02 
.76605+02 
.79705+02 
.79805+02- 
, .75975 + 02 
.8C005+02 
.80105+02 
. .80205 + 02 
.80305+02 
. .80405+02 
.80505+02 
-.80505 + 02 
’.80705 + 02 
^.80605+02 
. .80305 + 02 
. .81005 + 02 
.81105+02 
; .81205 + 02 
.81305+02 
. .ei405+02 
1.81505+02 
L. 81605+02 
.8 1 7 05+02 
» . 8l 8 0E+02 
‘.81905+02 
. .8200E+02 
'..82105+02 
.82205+02 
.82305+02 
.82405+02 
l .82505 + 02 
.82605+02 
' .82705+02 
.82805+02 
. .82905 + 02 
. .83505 + 02 
.831 05+02 
, .83205+02 
.83305+02 
.8340E+02 
I. 83505+02 
’ .83505+02 
: .83705+02 
.83805+02 
, .83955+02 
' .84005 + 02 
‘ .94105+02. • 
.84205+02 
I,. 84305+02, 
.84405+02 
.84505+02 
I .84 505+ 02 - 
,.84705 + 021: 
! . 84805+02 ■ 
.84905+02 
: .85005+02 
-.851 05+02 
* .85205+02 ; 
.85305+02 
.85405+02 
. .85505 + 02 
. .85605+02 , 
.85705+02- 
.85805+02 
’ .85905 + 02 
.86005+02 
• .86105+02 
.86205+02 
.86305+021 


H 


+3. 78552*33 
-0.74952+03 
-$.691 52 + 00 
-$.61 21 2 + 03 

-kttimri 

+3.451 3S+03 
0.4-1 292+33 
3.44532+00 
0.34352+33 

3 • 4 »2 $ 4 2+ 3 1 
0*46472+01 
f.'i-aBi £+oi 
3.485*2+32 
3. 45552+01 
3.40532+3? 
3. *51 3 72 + 3? 
0 .44792-01 
-$.*29522+00 
+0.oI342+33 

-:d. <55742+?: 

-$.*74232+0 3 
-$ .7 9 2 5 2 + C 0 
-$.*3215203 
-$.*33532+03 
-$.*34322 + 03 
-$.*33522 + 30 
-$.*32132 + 30 
-$.‘79752 + 33 
-$.'75152+00 
-$.71321+00 
-$.<53792 + 00 
-$•*33 732+30 
-$.4?:x72 + 30 
-$.£1552+00 
D.2905SO1 
3. ‘341 7200 
0.7203 23 0 0 

3.41412+01 
0*4 5 45 £ + 01 
0.43342+01 
w5v 13992+31 
0 .*16 7 32 + 01 
0.42212+31 
3.65972+00 
3.46392+00 

-$.*$922+0:: 

-$.45332+33 
-0.62 352+00 
-0.721 12+00 
-’0.73532 + 00 
-0.61 172+00 
-0. 63312+00 
-$.*3302+00 
-0.63712+00 
-0.32512+0? 
-$.<30542 + 33 
-0.77332+00 
-0.72532+00 
-0.66122+00 
-0.67012+00 
-0.41132+00 
-0. *27432 + 00 
-0.48522-01 
3. *24392 + 00 
5.60452+00 
3. *10172 + 51 
5.14 312+01 
5. *17532+01 
'-0. 49332 + 01 
0.47752+01 
5.43752+01 
5.62942+05 
0.30302+3? 
-0.41 322+03 
-0.45572+03 
-0.68 152+DD 
-0.69732+03 
-0.75522+03 
-0.63552+03 
-0. *3 2 9 32 + 53 
-0.63332+03 
-0. *33392 + 33 
-0.43322+03* 
-0. *31 2 32+ 0 3 
-0*78382+3? 
-0.71192+0? 
-0.63 2'3 2 + 0? 
-0.6 9 9 32 + 0 3 


-P.4550E+00 
-g. 49552+00 
-P. 52452+00 
-D. 64352+00 
-P. *54212+00 
-0.62612+00 
-P. 4 8902+00 
-D. 42592+00 
-P. 33142+00 
-P. '20382 + 00 
-P .32592-01 
0.4671E+00 
0.38142+00 
3.6809E+00 
3.73 37 E+ 00 
3.3013E+0*. 
0.7 9522+0 u 
0*721 52+00 
3.6315E+00 
0 .457 5E+00 
3.3255E+0C 
3.4 9552 + 00 
0.64 242-01 
-P.6 0 1 62-0 1 
-P.1 5212+00 
-0. *24182+30 
-P.31 952+00 
-D. *38582+00 
-P. 44392+00 
-P .48482*00 
-P.51 722+Orv 
-P.53722+00 
-P.51312+00 
-P.5325E+00 
-P.6021E+00 
-P.4171E+00 
-D. *36232+00 
-P.‘24 26E + 0 0 
-D .*35352-0 1 
0 .4 053E+00 
5. *31832+00 
• 3.5252E+00 
3.6937S+0C 
0.79U2+0C 

5 . *8 3 5 1 2 + 0 o 

0.74852+00 
0. 64132+0 1* 
0.5383E+O0 
3. *36732+00 
0.22922+00 
D. 99922-01 
-P.18 US-01 
-D. 42372+00 
-P.*2l 69E + 0 0 
-P.*2983E+00 
-P. *36752+00 
-D.4259E+00 
-D. 47312+00 
-P.5089E+OO 
-D.63252+0C 
-P. 54232+00 
-D. 53742+00 
-P. 51312+00 
-P.4558E+00 
-P. 39D22+0O 
-P. *281 12+03 
-D.4317E+G0 
5.4780E-01 
3. *25532+00 
•0.45842+00 

E:t?lSi:8S 

i-.mun 

3.67582+00 
3.54822+00 
3.4D31E+00 
0.'2584E+00 
3.43522+00 
3.1 473E-0 1 
-P. *91392-01 
-P.191 1E+0O 
-P.'2757E + 00 
-p. *34352+00 
-P.4131E+00 
-P.4605E+0C 
-P.4995E+00 
-P. 52592+00 
-D. 54132 + 00 


-0.07 962'+ 00 
-0. 4343203 
-0.47762+03 
-0.60912+03 
-0.6 27 62!+ 0 3 
-0.631 52+03 
-O.61R02+Q0 
-0.48372+03 
-0.42392+00 
-0.3340 2*+ 03 
-0. <2099 2*+ 03 
-0.61192-01 

§:3!”B8E 

0.52002+03 
0.66?7r+0? 
0.74172+03 
0.75092+03 
0.69962+30 
0.60502+0? 
0.48532+33 
0. *3551 2+33 
0.72462+0? 
0.4 O'US+D*' 
-0.45002-01 
-0.4 1902+03 

-0 .*2114 2+0 3 
-0.* 29712+03 
-0. <351 42 + 0? 
-0.41932+03 
-0.46602+03 
-0.601 02+03 
-0.62352+03 
-0.63182+03 
-0.62382+03 
-0.49592+00 
-0.44412+03 
-0. *35352+03 
-0. <24962+03 
-0.10082+03 
0. >79052-01 
0.27672+03 
0. 45892+03 
0.62672+03 
0. *72552+03 
0**75482+33 
D. 71962+03 
O.63572+03 
0.52172+3? 
0. *3 93 32+03 
0.76212+03 
0.43542+03 
0.4 741 2'— 0 1 
-0. 89892-01 
-0.43572+03 
-0.7 5972+03 
-0. *34232+03 
-0.43362+03 
-0.45352+03 
-0,49702+00 
-0.61822+0? 
-0.63092+00 
-0.62782+03 
-0.50612+03 
-0.461 82+03 
-0. <39032+03 
-0.78612+03 
-0.44742+00 
0.74342+01 
0.71882+00 
0.41532+03 
0.58602+00 
0.73372+03 
0.7532 g+00 
0.73542+03 
0.66372+03 
O.6567r+00 
0. 4311?+ 33 
0.79992+03 
C. 47132+00 
0.60542+01 
-0. 69982+01 
-0.45912+00 
-P.74662+00 
-0,32742+00 
-0.38692+30 
-0.44022+03 
-0.48702+03 
-0.61182+03 



- TABLE . 3 . CocfJTb .) 


i'.8b40E+02 -0 

f .8&50E + 02 -0 

1 .86605 + 02 -9 

> *8&7 0E+ 02 0 

v. 86802 + 02 0 

.8&90E+02 0 

:. 87002 + 02 D 

..87102+02 0 

'.8720E + D2 
V.8730E+02 0 

'.8740E + 02 0 

;.8750E + 02 3 

n .8760Et02 0 

.87702+02 -!? 

.8780E+02 -3 

.87902+02 -0 

.880OE+O2 -0 

.8gl0E+02 -9 

.8B20E+02 -'3 

'.8 8 3 0E + 02 -9 

.8840E + 02 -3 

'.8650E + 02 -9 

. .88602+02 -9 

‘.«S70E+02 -D 


•0.4471E+7D 

-9 
-9 
-0 
-9 
-0 
3 
3 
3 
0 
3 
*7 


• 4 8 4 9 E + 0 0 
.3*2902 + 03 
.'12312 + 3 ? 
.•15 25 2+ 0 0 
.49402+53 
.3947E+0* 
.•1317 2 + 01 
.16352+01 
.13932 + 01- 
.'13 3 92 + 01 
.1 535E+01 
.*93 9 92 + 0 0 
.44502+03 
.67312-92 
.'3 *2982 + 00 
.*5401 E + DO 
.*5 7 05E'+0D 
.*7 4 9 32 + 00 
.*79702 + 00 
.*3 2 392 + 0 0 
.*33702 + 00 
.*33972 + 00 
.33352+00 
.*31 34E + 00 


-P.5406E+00 
-P.'5 2?2E + 0C 
-P. 48202+00 
-P.'MSIE + OO 
-P.31 51E+00 
-D.1 B07E+00 
-P.-31 35E-02 
3.19452+00 
0.43862+00 
5. *50342 + 00 
0.7442E+00 
5.-3369E+00 
0.7890E+00 
0.*7075E+OC 
0.-5S69E+00 
. .3.44382+00 
0.*308fE+00 
0.1 7322*00 
3.48532-01 
-D.6421E-01 
-P.'l 6 45E+00 
-P .25 25E+00 
-P.-3237E + 00 
-P.3935E+00 


y,.8B90E+02 
;.8900E + 02 
,.891 0E+02 
.8 920E + 02 
. 8 9 30 E + 02 
'.8940E + 02 
...8950E + 02 
.89602+02 
'.6970E+02 
.8980E+02 
.8990E+02 
‘.9000E+02 


- 9 . 3 5 9 5 E + 3 3 
.•7554S + 0D -D.4 895E + 00 


.'701 52 + 3 3 
.62592+00 
.621 4E+0D 
.3*7 35E + 0 0 
.*13 5 IE + 0 0 
.6725S-P1 
.36902+03 
.*7 7 5 02 + 0 0 
.11972+31 
.*15912 + 01 
.4 3552+31 


-9.52332+90 
-P.6335E+00 
-P.5424E+00 
-P.5294E+00 
-P.4959E+00 
-D.4373E+00 
-P. -34812 + 00 
-P.2234E+00 
-P.61 242-01 
0.1 341 E+ 00 
0.3475E+00 


v 

-8 


0.6286 
0.6304 
0.6144 
0.4771 
0.4137 
0.3195 
0.4 909 
.281 1 
-.4 613 
0.3599 
0.6410 
0.6762 
0 .*7 463 
0.-7466 
0.6887 
0.690? 
0.4682 
0.337b 
0.2077 
0.8436 
-0.2923 
-0.4 317 
-0.2225 
-0.3017 

•0.4259 
•0.4711 
-0.6044 
-0.6251 
•0.5314 
•0.6209 
-0.4901 
■0.4347 
•0.34Q9 
•0.2 31 4 
• 0 . <7 3 0 5 
0.4 05 0 


E+DO 

e+oo 

e+oo 

s:+oo 

mi 

E+03 
2-01 
e:+oo 
2+0 0 
2+03 
E+oO 

e+qo 

e+oo 

E'+00 

e:+o? 
2+0 0 
e:+od 
2+00 
2-01 
2-01 
2!+0 0 
2+0 0 
2+0 0 

2+0 0 
2+00 
2+00 
2+0 0 
2+0 0 
2:+oo 
2 '+ 00 
2+0 0 
2+0 0 
2+0 0 
2-01 
2+0 3 


MP.-re: /yT-£*MEJ>,*T£ STEPS HAVB BE£M OMSrra*. 



-7716 L£ 4-, A -TRAJEtTOfLy cALCOLATEb FOR. 

Tva/O-T EmPERATURE Mop EL WITHOUT DELAYED rJEUTROfiJS 

4.530000 , 3. <33033030 + 03, 2.300000 „ 0.5300000 , -B 

TOL= 0.4 00-132 
INITIAL COOiriDSS: 

0.20000+ oo o. ioooe+oo o.ioooe+oo 


rcn 


5f (2) 


U3) 


.40000+ 0 0 
.2DOQ0+OO 
■*■ . 3 30 0 E + 00 
.43300+00 
* .500CE+0? 

.633DE+00 
. .7300E+0C 
.8000E+30 
.9D0CE+30 
, .1 0000+01 
.■•llOOE+01 
; .12000+31 
, *1 3005+01 
.1 4D0E+ 0 1 
*.1500E+01 
, .1&0OE+01 
^.1 7330+31 
*.1&00E+01 
.19030+91 

V.2000E+01 

-:iiS3i:»S 

;:3l*8USi 

. S«'2500Et 01 
i| . 2 6 0 0 0 + 0 1 
**.27000+01 
|'.2800E+01 
f .2903E+01 
,>30900 + 01 
1.31390+01 
> .32030+91 

i. 33 O'OE+0 1 

.'.340CE + D1 

.3500E+01 
.3600E+01 
*. 3 7 0 0 E + 0 1 
>38330+31 
^.39000+01 
3.4000E+01 
^.41 00 0+ 01 
S.420QE+91 
.43000*91 
.4400040! 
.45000*01 
i'*4600E*01 
i" . 4 7 0 0 E * 0 1 
: .4800 E + 01 
.49000+01 
■i .50300+0! 

I .51000+01 
£ .5200E+01 
*i .53000 + 01 
j .5 4 0 C E + 01 
1 .5500E + 01 
: .56000+01 
I .5709E+01 
! .580&E+01 
i .59000+01 

j. . 600 00+01 
!.6l00E+01 
-.62300+01 
, .63000+01 
; .6400E+01 
* .6500E+01 
f .6600E + 01 
■ s .670CE + 01 
, .68000 + 01 
4 . 6 9 3 0 E + 0 1 

.70000+01 

..710QE+D1 


0. '14 29S + 03 
0.35500-^1 
0. '29 3 30-01 
-’O.*2S4 2 0-Ti 
-0. '7 7350-0! 
-0 .42250+0 
-:o.'i5 3 90 + o:« 
-’0. *1353040? 
-‘0. -2027040: 
+3. <20940+3? 
-O. '23530 + 0? 

f Ilg+B: 

-0.43330+3: 
-0.97320-31 
-O. '53530-01 
3. 'H 7 200 2 
3. ‘57132-01 
3. *11250 + 03 
3. '16350 + 00 

5. *2 33 3E + 00 

0. ‘23550 + 00 
3. *25110400 

i:m\ u< 

0.46330+00 
0.43790+03 
3.79520-01 
0.4623 0-0 1 
-0.41520-01’ 
-0.95030-01 
-0.44230+00 
-0.47970+00 
-0. *20530 + 30 
+3.22720+3O 
-0.22 7 20*00 
-O. *22 150 + 0 0 
-0. *20510 + 00 
-0.47920+00 
-0.44330+30 
-0.98530-01 
-0.45130-0 1. 
0.41370-Dl 
0.7239 0 -D I 
0.43220+30 

BrillllloD 

3. '25370 + 00 
0. *27550*00 
0.25920+03 
0. *24190 + 00 
0.49710403 
0.43920+00 
3.73330-M 
3.651 5S-P2 
-0.67 3 50-0 1 
-0 .44 550 + 03 
-0.46390+00 
-0.20160+00 
-O. ‘22790 + 00 
-O .'24 27E + 0 0 
-0 .*24 59 0 + 00 
-0.-23790 + 03 
-0.21 380+03 
.'0,463 30 + 33 
-0.44390+00 
-0 .995 2 0-01 
-0.42350-31 
0 .*21 4 0 0-31 
3 .67350-31 
0.45530+33 


0.10960+00 
3.41300400 
3.43270400 
3,6 3 3 3E-01 
3.69750-Dl 
0.45910-01 
0.2133E-01 
-P.4309E-02 
-D. 29220-01 
-J0.5256E-01 
-D. 72950-01 
-P.6947E-01 
-P.'lOl 4E400 
-P.40930400 
-D. 109304-00 
-P. 43100400 
-D. 93330-01 
-P. 75150-01 

-D.*2959E-03 
3.-2925E-01 
D. 57730-01 
3.63340-01 
0.4 332 E+ 00 
0.4156E+00 
3.4 223E+0C 
D. 12330+00 
0.41300+00 
0.93380-01 
D. 71 460-01 
0.45850-01 
0.48220-01 
-P. 96590-02 
-D .36920-0 1 
-D. 61590-01 
-D. 83350-01 
-P.13310+00 
-D. *11230+00 
-P.1132E+00 
-P.4 1320 + 00 
+0 »•! 1 1 7E+00 
-P.96340-01 
-P. 79330-01 
-P. 55350-01 
-P. 26540-01 
0. *53570-02 
♦ 0. 67430-01 
0.88290-01 
0.95290-01 
0.4! 53E+00 
0.4295E+0C 
0.43410+00 
0. 42990+00 
0.4 1 74E+00 
0.97370-01 
0.7296E-01 
0.4 4 45E-01 
0.44150-01 
-P.4625E-01 
-P.45230-01 
-P.7145E-01 
-P.93B0E-01 
-P.4 1 1 3E+00 
-D. 12320+00 
-P. 42890+00 
-P.42790+00 
-P. 11 990+30 
-D .4 0 530 + 00 
-P. 63430-01 
-P. 55980-31 


0.4 3 3 5 09- D 3 
0.1 0730+03 
0.4 3690+03 
0.1 3D50+O3 
0. *6 8 3 3 0-01 
0.71! 20+01 
0.63040-31 
0 . -2 b 4 5 0- 3 1 
0.4 7620-32 
-0. 22680-01 
-O.456d0-O1 

_ n £. t * r\ * 


-0 


0.65030-01 
0.628/ 05-31 
-0.95330+01 
-0.4 3290+03 
-0.4 3500+03 
-0.4 31 30 + 03 
-0.91840-01 
-0.76810-01 
-0.56750-01 
-0.32580-01 
-0.5591S:-02 
0.22570-01 
0. 50010-01 
9.74750+01 
0.94930-01 
0.40900+03 
0.4 1630+03 
0.4 15406-03 
0.4 0750+33 
0.931 30-01 
0.73330-31 
0.49640-31 
0.235 7 0-01 
-0.33490-02 
-0.29690-01 
-0.543 4 0-01 
-0.75560-01 
-0.93030-01 
3.10570+03 
0.41280+03 
0.4 1410+03 
.43930+03 
97770-01 
0.83490-01 
0.57830-01 
0.33830+01 
9.96480+03 
0. 29920+01 
0.59690-01 
0.66150+01 
0.10720+03 
D. 12130+03 
0.12740+03 
0.12520+03 
0,1 1510+03 
0. 98060-01 
O.75470+O1 
0.48990-01 
0. 23290-01 
0.69590+02 
0.77260+01 
0.63250+01 
■0.65780-01 
0.43330+33 
0.1 1660+03 
‘0.4 234 0+03 
0.4 2370+03 
*0.4 1730+03 
*0.4 041 0+33 
*0.94360-01 


=8 


.633 



TABLE 4 -. ( co * Tb -) 


» .72D0E+01 
*.73&0E+01 
; .74ooe+ai 
J.7500E+01 
i .7600E+01 
l .7750E+01 

;.: f ? 88i : 8i 

■:: S188 ! i81 

; .g43OE + 0.1 

I . 8530E+01 

r 81 

. .8&30E+C1 
.8903 Et 01 
. .9C00S+D1 
.91 30E+ 01 
; .9233E+01 
; .933DE+31 
. .9133E+01 
, .9533E+31 
, r .9630E+01 
.-.9730E+01 
. .9333E+01 
' ,99DCEt01 
. . 1 0 30E+02 

. .131 OE + 02 
•; . 1 0 20E+02 
.1030E+02 
i .10fOE+02 

■ .1050E+02 
; .136&E+32 
! .10735+02 
f-.10§0E + 32 
L.1090E+02 

J . 1100E+32 
? .Ill 0E + 32 

' L .1 120E+32 
;.ji3oe+02 
r .i 1 40E+D2 
f . 1 1 50E + 02 
i . 1 1 60E+02 
* .1170E+02 
. . 11B0E+32 
. 1 1 90E+ 32 
.12D3E+02 
i .1210E+02 
.1 220E+02 
1.1 230E+D2 

K. 1240E+32 
4 .1 250E+02 
i .1260E + 02 
i'. 127 OE+02 

'.1239E+32 
. .1 290E + 02 
130DE + D2 
. 1 3 1 OE+02 
-.1 320E+32 
1333E + 32 
*. .1 340E + 02 
v . 1 3 5 0 E + 3 2 
'. 1 3&0E* 02 
it. 1 37 OEt 02 
.1 380E+02 
: .1390E + 02 
. .1430E+32 
, 5 .1410E + 02 
.1420E+02 
...1 430E+32 
. .14435+02 
.1450E+32 
.146JE+02 
. 1 4 7 0 E + 0 2 
r .1480E+D2 
t .1493E+02 
. 1 5 0 0 E + 3 2 
. .1510E+02 

■ 1 520S + 02 
* .1 533E+32 
•.*.154 3 5+ 02 

.15505+32 

■ .1553E+32 
. .15705+32 

.15805+02 

4.15905+02 


0. '21 295 + 00 
3. '2511 5 + 03 
0. '29295 + 03 
0.303 9 s:+oo 
0.. 29225 + 93 
0.'25S5E + &0 
3. •20595 + 33 
3.44335+03 
3. '571 45*31 

: 3 : 88ii | 3 ? 

-:3 .'1-3 5 3 * + 00 
-0.49545+00 

-Oli?355+SS 

•0 .*25 5 + 

- Ov 2 ' 329 E + OC 
-3 .' 1 9 9 4 £> 3 3 
-* 3.1 5 5 1 £>03 
- 3 . 13 33 £>30 
-; 3.*3 7 7 4 E -:01 
D .' 3175 E -:31 
o.i 344 e*od 
3 . 1 753 E + D 3 
3 **24 3 3 EO 3 
3 .' 2 * 913 E >30 
3. ’3244 Et -33 

3.’3*33l£f33 
0 .* 3i59E : O3 * 
0 *'27 Si £>0 0 
3 • *21 51 £ + 33 
D.'U35E + 33 
3.5955E-D1 
>3.'2045E-Dl 
-0**54 36 £-01 
.’0 .45 3 3E> 0 3 
-P. *21 3 3 £>5 3 
-3.24 3 7 £+3 3 
-:0.'27UE + 00 
-*5.'2354E + 03 
-0.*2!353E+O3 
•0 **27 31 £>0 0 
-:0.‘2*430£+3Q 
.’0**21 0 9 £ + 0 0 
-P .1 6 23E + 0 3 
-0.4023E+D3 
-' O . WSE-fcl 
3*4233E-D 1 
0.1*21 5E + 03 
3.1994E+03 
3*'2595E + 03 
0**32 4 7 £ + 0 0 
0 .*36 B 4E>0 3 
0**36 5 4 £>0 0 
3.*3437E + 30 
3.*295l£ + 03 
0 .‘2253E + 00 
3.1>1 3£>33 
0.6179E-P1 
.:0.»3515£-X>1 
.:0.‘i*l40£>00 
•0 **131 0 £> O 0 
-0.'2344£>P3 
-!3.*27 35£>03 
-;0.*2935£t00 
-9 .<30 9 7£> 0 0 
.;j.*3 07 5£t00 
• 0 .'29 2 4 £>00 
•0 .'2*64 4£>30 
-p.'2'2 3 8 £>0 0 
p • *17 0 9 £> D 0 
0 .4 *3 5 2 E>30 
0. *31 2 5 £-91 
0.61 5 3 £-91 
3 • 1 3 3 4 E + 0 3 
0**22 3-5 S>00' 
0. *30 0 0 £>00 
3.<3597E + 03 
3.‘3949E>33 
0 *>4 D 3 l £ + 33 
3 . *3 7 3 3 E + 0 3 
0.*31S2£ + 03 
0 . 2 3 S 2 £ ♦ 0 0 
0.142 3 £>03 
0.441 0£-9l 
-0 . *4 9 9 9 £ -9 1 


«D.*2396E-01 
0.13B9E-01 
0.4632E-01 
0.7974E-01 
0.1035E+00 
0.43030*00 
3.1433E+00 
3.145/0+00 
3.1 4O3E+0O 
0.1 2510+00 

knmm 

0.4230E-01 
3.97100-02 
-D.'i3l OE-Ol 
-l0.o4 05E.0l 
-D .3 1 7 9E-0 1 
-P .4 C 5 20+0*3 
-B.1 2320+00 
-0.1 5 5 2 E > 0 0 
-9.1 403E+00 
-D.1392E+0C 
-D.1 ?3 7EtO 0 
-D . 1 1 1 5E+00 
-0 .*3 7 3 1 E* Cl 
-D.5675E.01 
-0.2129E.O1 
0.173BE.C1 
0.6532E.01 

0.-9202E-01 
0.1227E+0O 
3.U54E+00 
3.1581E+00 
0.16D3E+00 
3 . 1 5 1 3 E 1 0 0 
0 .1 332EtO 0 
3.1073E+00 
3.75B3E.01 
5.4105E-01 
0.5055E.02 
-D.-3026E.01 
-D.6327E.01 
-D.9258E-01 
-D.1170E+00 
-D.1357E+DC 
-D.1477E+00 
-D . 1 5 ? 4 E 1 0 0 
-0.1493000 
-0.1 33iE«D0 
-0.1 1 38E+00 
-D.9175E-01 
-D.573BE.0l 
-0 .1 B75E-0 I 
0.*2342E-01 
3.6575E-01 
0.1050E+00 
0.1378E+00 
0.1S15E+0O 
0.1 7 39000 
3.1743Et00 
0.1632E+00 
0.1419E+00 
0 .1 1 25E>00 
0.7 75 4E-0 3 
3.«3943E-01 
D.6895E-03 
-0.5755E-O1 
-D .727 1 E-D 1 
.0*1 93 7 £>00 
-0 .1 294E+00 
-P.1437E+00 
-D.1609E’f00 
-0.*1 653Et3O 
-0*1 61 2E+00 
-0.1435E+OO 
-D.1259E+09 
-0.97O1E-91 
-0*5959E-01 
-0 *1 6 S6E“0 1 
3*‘29SOE-01 
3. 7587E-91 
0.1 135E+00 
3*1538 Et 0 0 
0 .1 735E+0O 
0.1 938E+00' . 
0 .1 8 97 El 00 
3*17 6 0 E> 9 0 
0**1 51 4E>00 


- 0.5 89 3 £*-01 
• 0. *2889 £-01 
0 * 4081 2*02 
0 . 6791 £-01 
0**701 9 £'» 8 1 
0 . 9846 E -01 
0 . 1 2 0 5 £’-t 0 3 
0.1 34 5 £>00 
0.1 3 9 7 £> 0 0 
0 . 1 35 6 £>03 
0 . 1 2 3 1 £> 0 3 
0.1 0 3 2 £> 0 3 
0 . 7758 5-01 
0 . 481 4 £-01 
Q .1 67 0 £’- 0 l 
- 0 . 1 498 E -01 
- 0.4531 £-01 
- 0 . 7 2 9 3 £» 0 1 

- 0 . 9 6 5 5 £- 0 1 
- 0 . 1 152£>03 
- 0.1 2 3 1 £> 3 3 
- 0.1 3 4 5 £> 3 3 
- 0.1 3 3 9 £> 3 3 
- 0.1 2 6 3 £’> 3 3 
- C.l 1 0 8 £> 0 3 
- 0 .^ 360 £-31 
- 0 . 631 9 £-31 

- 0 . -2 5 ft 9 £ - 0 1 
0.*9449£'— 02 


0.4645 
0.6143 
0 . 1 117 
0.1 347 
0.1 48 7 
0.1 577 
0.1467 
0.1315 
0.1 385 
0.7979 
0.4725 
0.1294 
- O .0128 
- 0.5373 
- 0.6300 
- 0.1 078 
- 0.1 272 
- 0.1 403 
- 0.1 464 
- 0.1 449 
- 0.1 355 
- 0.1 182 
- 0.9336 
- 0.6182 
- 0.0503 
0 .1 49 1 ; 
0.5536 
0 . 9329 : 
0.1 257 
0.1 498 : 
0.1 63 7 
0 . 1666 : 
0.1585 
0*1405 
0.1144 
0.6233 
0.4651 
0 . 9201 ! 
- 0.0766 
- 0.6235 
- 0.9338 
- 0.1 195 
- 0.4397 
- 0*1531 
- 0.4589 
- 0.4565 
- 0.4458 
- 0 * 1 264 
- 0.9391 
- 0.6411 
- 0.<236 4 
0.201 8 
0.6438 
0.4056 
0.4 104 
0.4658 
0.1798 
0.4 815 
0.4712 


£-31 
£-01 
£>0 3 
£> 0 D 
£>0 3 
£> 0 3 
£>0 3 
£>00 
£> 0 0 
£-01 
£-01 
£-01 
£-01 
£-01 
£-01 
£>03 
£>0 0 
£>00 
£>0 0 
£>00 
£>00 
£+0 3 
£-01 
£—01 
£>01 
£-01 
£-01 
£-01 
£>00 
£> 0 0 
£>00 
£>00 
£>0 0 
£>00 
£>00 
£-01 
£>01 
£>02 
2>01 
£-01 
£-01 
£>0 0 
£>0 0 
£>0 0 
£> 0 0 
£>0 0 
£>00 
£>0 0 
£-01 
£-01 
£-01 
£>01 
£-01 
£>00 
£>0 3 
S>0 3 
£>0 3 
£>03 
£>0 3 



TABLE 4-, (coA/TD.) 


:•.} 045E + 03 
...1046E + 03 
t.}0*7E+03 
^.104BE*D3 
ji.4 049E+03 
f;.1050E*D3 
S..1051E+03 
' .10525+03 
;-.10§3E + 03 
.;.lD54E + 03 
.VI 055E+03 
tf.lQ56E + D3 
. .1057E+03 
V.1058E+03 
>.10535+03 
r.lOSOE+03 
.10615+03 
...1052E+33 
■..1063E+03 
*.1064e+(J3 
’.10555 + 03 
. 1 0 6 6 E + ? 3 
.1067E+03 
i . 1 0 6 8 E + 0 3 

’.1071E+D3 
V1072E+33 
r.i 07 3E+ 03 
’.1074E + 03 
V 1 07 5 E + 03 
^.1076E+03 
-,■»*! 07 7 E + 03 
..1078E + O3 

,4.1 079E + O3 

VI 080E+03 
5.108 1 E + 03 


,.1085E+03 
. .1086E+O3 
-.10875 + 03 
*.}08BE+03 
..1089E+33 
4-.1090E + O3 
.1091E+33 
-.10925+33 
.1093E+03 
.1094E+03 
.t . 1 0 9 5 E + 0 3 
, .1096E+03 
W.1097E+03 
i-. 109 8 E+03 
K. 1099 E+03 
i-.‘l 1 0DE + 03 

I..U01E+03 
_ .1 1 02E+O3 
» .1 1 03 E+03 
UU?4E+03 
y .1 1 &5E+33 

t:im§;8! 
t:l 1S8USI 

l .HtOE + 03 
t.llIlE+03 
V.1112E+03 
V1113E + 03 

V. 1114E+03 
s .1115E+03 
fc.l 1 1 6E+03 
. .1117E + 03 
i, . 1 1 1 8 E + 0 3 
„.1119E+03 
i .1 120E+33 
E.U21E + 03 

W. 1122E+03 
*.ll23E+03 
■■ .1 1 24E + 33 
» .1125E+03 
. .11 26 E+03 
VI 1 27E+03 
S.1128E+03 


-0.4539E+OD 
»S> .'30 365 + 00 
-;3.fD37E-01' 
0 • 1 8 5 6 E + 3 0 
3. ? 3 3 35+0 3 
0.9275E+0D- 
0 . 1 -3 3 9 E + 0 1 
3.'1535E + 01 
wO.'IS 5 lE+Dl 
O.-1773E+01 
0. '14335 + 01 
3. 9 1 7 45+33 
3. '391 45 + 03 
-’3 .'4 3 925-Di 
-0. ’34315 + 33 
-3. 64 3 15+33 
-’3.67235 + 03 
-’3. '74 3 75 + 33 
-’3. '79 4 35 + 33 
-3.61 975+33 
-;3. 631 75 + 33 
-.'3.63345+33 
-0.62595+33 
-:3 .63925+ 33 
-’3. ’731 95+03 
-:3.'7 41 3 5 + 33 
-’3.63375+33 
-:3. 63335 + 03 
-0. '49255 + 0? 
-!3.'34i 75*33 
-0 .'1-3 995 + 03 
0. *1-2365 + 33 
0. '45425 + 30 
0.64395+03 
0. '1 2595 + 31 
3. <16295 + 31 

-kmmn 

0 . 15252*01 

0.10322+01 
D.49542+00 
0. *37 5 02 -01 
-■3. ’292 62+00 
542+00 
+0. *55 302+0 0 
-0.73712*00 
-0.7375E+00 
-0 .‘3t 552+0 D 

-0.0311 s:*oo 

+0. *93472*00 
-0.029 22+00 
-0 .*3 1 4 5 £ *03 
-p . 73 972*00 
-*0.75232*00 
-*0 *69 3 92 + 0 3 
-10.62*12*03 

zkmum 

+3. 10972*03 
D .59352"Pi 

0 .'37 532+30 
3. ‘7557203 
0.11712+01 
0.15502+ 31 
0.10242*01 
19592 + 01 
0.16152+01 
3*4 1512+01 
-'3*6-1 312 + 30 
0**1*27124-03 
-13.7'3432 + QO 
-’0.4 7 5 52+ DO 
-0 *6291£>0D 
-10.72312 + 00 
-0**77952 + 03 
-p. 01 252+03 
-13 .6 29 72*3 3 
-'3*6355 2* 3 3 
-13.6 32 3 2,4-0 3 
-p.*3 1 9 42+3 3 
-!3. ‘79712 + 03 
-*3. ‘75292 + 03 
-0 • *7 1 3 5 2 + 0 3 
-0.64422+33 
-0.64312*03 
-0*41542+33 
-0. ‘23372 + 03 
-0 .'34 9 32 "0 2 


-P.5322E+00 
-0.51 1 82+00 
-D. 46932+00 
-P.09 98E+00 
-D**298 2E + 00 
-0**16352+00 
0.1 2982-01- 
0.71 42E+00 
0.42342+00- 
0.6097E+00 
3**74022+00 
0*7 9 342+ OG 
0.7&9OE+OO 
3*68452+00 
3.65372+00 
3.4272E+00 
0 .‘28 3 8 2+ OG 
0*15552+00 
D. <34232-01 
-D. 76102-01 
-D. 17422+00 
-D. *26322+00 
-D. *33452+00 
-0.3 97 4 E+ 00 
-P. 44912+00 
-0.48962+00 
-P.51 3 32+00 
-P. 53412+00 
-D .5 3 5 2E*00 
-D.51 3BE+00 
-p. 48122+00 
-0.41 762+00 
-P.'3?29E+00 
-D. 19272+00 
-P. <25842-01 
3. *17142 + 00 
3. <391 BE + 00 
3.67542+00 
3.72132+00 
0.79142+00 
0.7819E+0O 
0.7076E+OO 
3.6920E+00 
0.4570E+00 
0*017 9E + 00 
3.48352+00 
3.58912-01 
-0.54142-02 
-P.l 5 4 92 + 00 
-0.24352+00 
-D.62 3 3E+ 0 j 
-0.00562+03 
-P.4 59 8E+00 
-0.43272+00 
-0.51 402+00 
-0.53282+00 
-P.53732+00 
-P.52502+00 
-P. 49242+00 
-D. 43492+00 

-P.64732+00 
-D.-2249E + 00 
-P.65522-01 
3.12552+00 
3.33562+00 
D.5332E+00 
3.69712+00 
0. 78502+90 
5.79222+00 
D.7299E+00 
0.621 OE+O0 
3.4683E+00 
3**34382+00 
0. <21282+00 
0. *35552-01 
-P.-3029E-01 
-0.13392+00 
-D. *22522 + 00 
-P.*30 47E+0 0 
-P.0725E+GO 
-P.42922+00 
-P.4747E+00 
-P.5337E+00 
-P. 5335E+DO 
-P.53942+00 
-D .53322 + 00 
-0.63262+00 
-p.451 3E+0C 


-0*6 22 3 £+03 
-0.6221 2*03 
-0.60432*03 
-0.46472*03 
-0.39892*03 
-0.3 0222* 0 3 
-0.17162*03 
-0.916025-02 
0.17Q82*03 
0.37382*03 
0.64762*03 
3.67402*03 
0.73542*03 
0.72972*03 
0.66832*03 
0.66852*03 
0.44742*03 
0.31842+33 
0.19062*03 
0.69582-01 
-0.4)772-01 
-0.1 4232*03 
-0.73082*03 
-0.3 3802+03 
-0.37392*03 
-0.42842*03 
-0.47162+03 
-0.63282*03 
-0.52112+03 
-0.62462+03 
-0.61072*03 
-0.47602+03 
-0.41592+03 
-0.32582*03 
-0.7021 2+03 
-0.44632-01 
0.1 3992*00 
0.33532+03 
0.61622+03 
0.65472*03 
0.73042*00 
0.73802+03 
0.69652*03 
0.69292*33 
0.47462*03 
0. 34612+03 
0.71732*03 
0.94272+01 
-0.1 9422*01 
-0.12222*00 
-0. *21352 + 00 
-0.79322*03 
-0.36152*33 
-0.41852*03 
-0.46422*03 
-0.49812*03 
-0.61932*00 
-0.62622*00 
-0.51642*00 
-0.48652*00 

-0.43232*03 
-0.34932*03 
-0.73272*03 
-0.62072+01 
0.97972*01 
0.79342*00 
0.49042*00 
0.63082*03 
0.72152*03 
0.74372*00 
0.70372*00 
0.61742*03 
0.53292*03 
0.37532*03 
0.74572*03 
0.12092*00 
0.47762+02 
-0.13062*03 
-0.1 9452*03 
-0.77662*03 
-0.34772*03 
-0 .4 07 32*0 D 
-0.45562*03 
-0.49222*03 
-0.61642*03 
-0.62682*03 
-0.621 32+03 
-0.49612*03 



TABLE A. (cofi/Tb.) 


£f*l 129E+03 
$Vll30E* 03 

KllliltSI 

i:Jlll|}§3 

HrlllJRj! 

MM5E+03 
!r.H46E*03 
.".'ll 47E+03 
h p.U4 8E + 03 
+.1 1 49E*33 
: .1 1 5OEt03 
i.ll 51E+33 
. .1152E+03 
v'.ll 53E+03 
.1154E+03 
..1155E+03 
, .1156E+03 

*.ll59E+03 

i.ll50E+03 

frrti tieJ81 

g-.Jl&3E*03 
1.U64E+03 
fe.'l 1 65E+03 
Ivi 1 66E+ 03 
lMl67E*03 
*>.4 1 68E + 03 
j* • ! 1 69E+03 
. .117PE + 03 
.11712+03 
.1 172E+33 
.4173E+03 
.11 74E+D3 
. .1 1 7 5 E + 03 
V.1176E+03 
11 77 E+03 
.-.1 178 E+03 
".*1 1 79E + 03 
. 1 1 8 0 E + 0 3 
-.'ll 51 E+03 
.U82E+03 
. .1 183E+03 
;.1lB4E+03 
: .11S5E+03 
.1 186E + 03 
: .11 87 E + 03 
, 1 1 8 8 E + 3 3 
-.1189E+03 
1 90E+03 
... 1 1 9 1 S + 0 3 
; .1192E+03 
,‘.11 93 E+03 
I • . 1 1 9 4 E + 0 3 
r; .l 195E + 03 
.1 196E + 03 
.1 1 97E+03 
. .11 98E+03 
.1199E+03 
.1200E+03 


0 .'29 782+33 
3.-S543E + 33 
0 • ‘1 0 7 7 E + 0 1 
3.-14912 + 31 



0. '12 7 3 £ + 01 
0.'/ 3 39E + 00 
3. 77 7 8 E+03 
-;3.'lo36S + 03 
-0 •'|311 St 0 3 
-iO.-dOIdt + OO 



-!3 .'3 3 7 2. 

•;3 .>3 2 7 5E + 0 3 
-0 • -3 3 5 7 E + 3 0 
-.'3.-3 3 42E + 33 
-p . -3 2 3 3 E + 3 3 
-V .- 3 3 39* +33 
-0 .'7 7 23 E+3C- 
-10.72 7 42 + 3D 
-;3 . '5 5 3 5 S + 3 3 
-0 .'5 7 4 4 2*33 
-’3 . 4 5 2 1 E + 3 3 
-0 ,'2’3 5 3 2 + 33 
-!> • *85 5 4 — P 1 
3.-21 99E + 03 
3.-57 27 E+0 3 
3.-9 7 3 2 E + 03 
3.-13922+31 
3.1 7 3-3 E+ 01 
v-3 .4-8 35E + 0 1 
3.17732+31 
3.13932+31 
3.-35 49E + 5 3 
3. -3-3 95 S + 3 0 
-p. <3 2 9 32-Pl 
-13. '3 -7 3 3 E + 03 
-p. -55332 + 33 
-;:y.63 572+33 
-;3 . -7 5 3 3 E ♦ 0 D 
-p. -33332 + 33 
-'3.-3245 2+ 3 3 
-P.-3 352E + 33 
— p »*3 3 d 3 E + 3 3 
-P .-32 7 7 2+3 3 
-0.-31 322+33 
-13.-7 3 21 E.+ C3 
-’3 .-7 4 3 SE + 3 D 
-;3 .-S3 1 SS + O 3 
-p.‘5 9 9 4E+0 3 
-0.-48 51 2 + 0 3 
-13 .-3 -3 2 3 E + 03 
-0 .1 -2 5 5S+ 03 
0.4 4352*30 
3 .4 3 3 5E ,+ 0 .3 
3.-37552 + 03 

5 .1 -2 9 52'+ 01 

3.16522+31 
^.‘18722 + 31 
3.-18 2? E + 31 
3. 1-3 3 72 + 01 
0.4001 2 + 01 
0.451 3E+ 03 
D .*32 3 4.5-3.2 
-'3.-31 7 5E + 30 
-;3 . -5 3 3 7 S + 0 3 
-;3 . -5 6 3 5 E + 0 3 
— 13 .*7 4 4 55 + 00 
-]3 . -7 9 2 9 2+ 3 3 


-P.3710E+00 

-0.255BE+O0 

]8 86E + 00 


. _955e + 00 
0.6675E+00 
0.-/736E + 00 
3 .799OE+O0 
0.7506E+00 
0.5499E+OO 
3. -5 207 El + 00 
0.-381 4E+0O 
0.2437E+00 
O.jl 40E+OC 
-P.465 /E-02 
*0.11 12E + 0 0 
-0.-2O54E+00 
-0.-28 7 6E+&0 
-P.-3532E + CC' 
*0 .41 7 5E + 0 0 
*■0 . 4 5 5 5 R + 0 0 
-0.5O24E+00 
-P .-52 7 2E + O0 
-0.-5335E + O0 
-P.5344E+G0 
-0.51,1 8E + C0 
-P.4665E+C0 
-0.-3936E + OO 
-0.-2877E + OO 
-D.-1452E + 00 
0. -3 323 E - 01 
0.-23 8 3 E+ 0 0 
D.448BE+0O 
3.6327E+00 
O.-7558E+O0 
0.8015E+00 
0.-7690E + 00 
0.6783E+OO 
O.5537E+O0 
0.41 53R+0O 
D.-2753E+0O 
O . 1 4 4 2 E + 0 0 
D.-2 2S3E-01 
-D.-36 9 3E-U1 
-0.18 41 E + 00 
-P.-2692E + 0O 
-D.-3426E+00 
-0.4O45E+DO 
-D.4555E+0U 
-P.4951 E+00 
-D.5228E+00 
-P.5376E+O0 
-D.5376E+00 
-P.-5198E+00 
-P. -48352 + 00 
-P.41 48E+0O 
-P.-31 7 3 E+0 0 
-P.’l 839E + 00 
-P.-135 9E-01 
3. 1357 E + 00 
3.-3989E+O0 
3.5928E+00 
3.-7344R+0O 
3.-7993E + 0C 
0 .7 8 43E + 0 0 
0 .70545+00 
O.-58&9E+0O 
3. 4 532 £+00 
3. *31 335*00 
3. 1 759E+O0 
3.-5I 47E-01 


-0.4478 EH- 00 
-0.-371 6 E+00 
- 0.262 9 g:+o 5 
-0.1 198iH03 
0.546921-01 
0.C488S+03 
0.44062+03 
0.60222+00 
0.- 70812+00 
0.7460 2+ 03 
0.-71892+00 

8:IHfp83 

8:8939388 

O.-l 491 2+0 0 
0.-30682-01 
-0.7? 372+01 
-0.1 7402+00 
-0.75912+00 
-0.-33262+00 
-D. -39492+03 
-0.44582+03 
-0.4 85 3 E+00 
-0.5125 E+00 
-0.5 26 3 E+03 
-0.524 7 E+03 
-0.5346 EH 03 
-0.4623E+0D 
- 0 . -3 9 3 1 E+ 0 3 
-0.-2923E+0D 
-0.1 5722+03 
0.4076E:-D1 
0.70212+03 
0.39732+03 
0.56892+00 
0.69012+03 
0.74452+03 
0.73182+03 
0.66452+03 
0.55062+00 
0.4371 EH 03 
0.-30702+0 3 
0.1 789E+03 
0.581 5E:-D1 
-0.52632-01 
-0.15212+03 
-0.74002+00 
-0.-31632+00 
-0.-38132+03 
-0. 43532+03 
-0. 47732+03 
-0.50762+03 
-0.52492+00 
-0.52722+03 
-0.51232+03 
-0.47552+03 
-0.41322+03 
-0. -32052+03 
-0.1 9382+03 
-0.-33182,-01 
0.-1541 2+03 
0.-351 02+03 
0.63142+33 
0.66722+03 
0.73862+03 
0.741 62+00 
0.58602+33 
3.58922+33 
0.45932+03 
0.-33932+03 
3.-21012+03 


No-rE ’ 


ibJTE.R.MEblAT£ STBPS HAVE. BBBfV OMt~T"T££> 4 
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